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I-— i! ABSTRACT 

Aims. Our current knowledge of high-mass star formation is mainly based on follow-up studies of bright sources found by IRAS, and 
Q , is thus biased against its earliest phases, inconspicuous at infrared wavelengths. We therefore started searching, in an unbiased way 

5^ . and in the closest high-mass star-forming complexes, for the high-mass analogs of low-mass pre-stellar cores and class protostars. 

' Methods. We have made an extensive 1.2 mm continuum mosaicing study of the Cygnus X molecular cloud complex using the 

, MAMBO cameras at the IRAM 30 m telescope. The ~ 3°^ imaged areas cover all the high-column density (Ay > 15 mag) clouds of 

' ' . this nearby (~ 1.7 kpc) cloud complex actively forming OB stars. We then compared our millimeter maps with mid-infrared images, 

' and have made SiO(2-l) follow-up observations of the best candidate progenitors of high-mass stars. 

" ' ' Results. Our complete study of Cygnus X with ~ 0.09 pc resolution provides, for the first time, an unbiased census of massive 

^ , young stellar objects. We discover 129 massive dense cores (FWHM size ~ 0.1 pc, Mj 2 mm = 4 - 950 Mq, volume-averaged density 

\l . ~ 10' cm"''), among which ~ 42 are probable precursors of high-mass stars. A large fraction of the Cygnus X dense cores (2/3 of 

the sample) remain undetected by the MSX satellite, regardless of the mass range considered. Among the most massive (> 40 Mq) 
cores, infrared-quiet objects are driving powerful outflows traced by SiO emission. Our study qualifies 17 cores as good candidates 
for hosting massive infrared-quiet protostars, while up to 25 cores potentially host high-luminosity infrared protostars. We fail to 
QQ ■ discover in the high-mass analogs of pre-stellar dense cores (~ 0.1 pc, > lO'* cm"^) in Cygnus X, but find several massive starless 

' clumps (~ 0.8 pc, 7 x 10^^ cm"^) that might be gravitationally bound. 

' Conclusions. Since our sample is derived from a single molecular complex and covers every embedded phase of high-mass star 

, formation, it gives the first statistical estimates of their lifetime. In contrast to what is found for low-mass class and class I phases, 

• • . the infrared-quiet protostellar phase of high-mass stars may last as long as their better-known high-luminosity infrared phase. The 

' statistical lifetimes of high-mass protostars and pre-stellar cores (~ 3 X 10* yr and < 10^ yr) in Cy gnus X are one and two order(s) of 

' magnitude smaller, respectively, than what is found in nearby, low-mass star-forming regions. We therefore propose that high-mass 

r> I pre-stellar and protostellar cores are in a highly dynamic state, as expected in a molecular cloud where turbulent processes dominate. 

Key words, dust — H n regions — ISM: individual (Cygnus X) — ISM: structure — stars: formation — submillimeter 



1. Introduction 

High-mass (OB, > 8 Mq) stars, though few in num- 
ber, play a major role in the energy budget of galax- 
ies. Our current understanding of their formation, however, 
remains very schematic, especially concerning the earliest 
phases of the process. High-mass stars are known to form in 
dense cores within molecular cloud complexes, by accretion 
jYorke & Son nhalter 200?; 'Kru mholz et al. 20071 1 and/or coa- 
lescence (e.g. Bonnell et al. 200 ft . The copious UV flux emit- 
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ted by a newly-formed central star heats and ionizes its parental 
molecular cloud, leading to the formation and development of a 
hot core (e.g. IHelmich & van Dis hoek 1997) and afterwards an 
H II region (see a review bv lChurchwell 19991 1. 

In the past few years, the progenitors of ultracompact H 11 
(UCHii) regions (sources that have not yet begun to ionize 
their environment in a detectable way) have been searched 
among IRAS point sources. The criteria generally used aim at 
selecting high-luminosity (> lO-' L©) stellar embryos embed- 
ded in a massive envelope and associated with hot gas but 
no UCH II region (see IKurtz et al. 20001 1. In fact, most sur- 
veys focussed on red IRAS sources with colors satisfying the 
[Wood & Churchwell (1989)| criteria (originally set for selecting 
UCH II regions) but not detected in radio centimeter surveys. The 
presence of an envelope around the IRAS sources is confirmed by 
their association with high-density gas, and the existence of hot 
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gas by the detection of a hot core and water or methanol masers 
(e.g. lBronfman et al. 19 96 ; Plum e et al. 19971 1. In this way, stud- 
ies by e.g., IMolinari et al. ( 2000)1 and [Sridharan et al. (2002)| 
have identified massive young stellar objects harboring high- 
luminosity infrared protostars (see also IMueller et al. 20021 
IFaundez et al. 2004l i. 

Such surveys are biased against younger and, thus prob- 
ably colder, massive young stellar objects which could be 
inconspicuous in the mid-infrared bands of IRAS and the 
Midcourse Space Experiment (MSX). These high-mass analogs 
of low-mass class protostars and pre-stellar cores (cf. 
lAndre et al. 20001 1 should be best detected via far-infrared to 
(sub)millimeter dust continuum and high-density molecular 
line tracers. Serendipitous discoveries of a few infrared-quiet 
massive young stellar objects have been made in submillimeter 
continuum maps (e.g. Molina ri et al. 199^ IMotte et al. 20031 
ISandell & Sievers 20041 IHill et al. 20051 1 and as absorp- 
tion regions in mid-infrared images (often called "infrared 
dark clouds", e.g. ICarey et al. 2000} IThompson et al. 2005 



IRathborne et al. 20061 1. These studies suggest that, in the high 



mass star formation process also, infrared-quiet phases must 
exist and be associated with significant populations of young 
protostars and/or pre-stellar cores. The lack of statistics, homo- 
geneity, and angular resolution of studies such as the above, 
however, prevents the determination of the basic characteristics 
of these infrared-quiet objects. 

To make significant progress, one needs to search in a sys- 
tematic way for the earliest phases of high-mass star forma- 
tion in nearby complexes. One of the most appropriate meth- 
ods is to survey entire molecular cloud complexes in which 
high-mass stars are forming, using tracers of high-density clouds 
and protostellar activity signatures. We therefore started a mul- 
titracer study of Cygnus X, which is one of the richest molecu- 
lar and H ii complexes located at less than 3 kpc from the Sun. 
The molecular cloud complex is massive (4 x 10^ M©) and ex- 
tends over ~ 100 pc in diameter (cf. jLeung & Thaddeus 19921 ). 
Found to be part of the local spiral arm, this complex has 
been located at 1.7 kpc from the Sun by the recent '^CO 
study of Schneider et al. (20 06) Cygnus X has a rich collec- 
tion of H n regions (see centimeter free-free emission maps 
by IWendker et al. 199T1 and infrared images from MSX pre- 
sented in Figs. 7-9 of Schneider et al. 2006), indicating re- 
cent high-mass star formation. It is associated with several OB 
associations (see ,Uyaniker et al. 2001 and references therein), 
amongst them one of the largest in our Galaxy (Cyg OB2, cf. 
IKnodlseder 20001) . It also contains several well-known massive 
young stellar objects like DR2 1 , DR2 1 (OH), AFGL 259 1 , S 1 06- 
IR, and W75N (e.g. Smith et al. 1998; v an der Tak et al. 19991 



[Schneider et al. 2 002: Sheph erd et al. 2004| . Small parts of the 
Cygnus X cloud complex have already been observed in 
(sub)millimeter continuum emission (e. g. [Chandler et al. 1993al 



[Richer et al. 1993l|Vallee & Fiege 20061 IDavis et al. 2007l l. 

The present paper tackles the following key questions: do 
high-mass pre-stellar cores exist? What are the lifetimes of 
high-mass pre-stellar cores, infrared-quiet protostars, and high- 
luminosity infrared protostars? What are the main physical pro- 
cesses leading from pre-stellar objects to high-mass protostars? 
We report an extensive 1 .2 mm continuum study of the Cygnus X 
complex complemented by SiO(2-l) follow-up observations of 
the best candidate progenitors of high-mass stars. From the 
MAMBO-2 imaging (~ 3°') of the entire complex presented 
in Sect. |2l we make a complete census of the compact cloud 
fragments and larger-scale structures (Sect.l3]l. We determine the 
main characteristics of these new millimeter sources in Sect. |3l 




Fig. 1. Extinction map of the Cygnus X complex derived by 
Bontemps et al. (in prep.) from the stellar reddening of back- 
ground stars in JHK, using the 2MASS database and a pixel size 
of 1.3'. The high-extinction regions are yellow and outlined by 
the 15 mag contour level. The three molecular cloud ensembles 
identified by Schneider et al. (2006) are labeled and the nearby 
Cygnus OB2 association is indicated. The fields mapped with 
MAMBO or MAMBO-2 are schematically outUned and add up 
to - 3°\ 

and search for signposts of protostellar activity including SiO 
emission (a tracer of outflow activity) in Sect.U Section|5]iden- 
tifies 129 massive dense cores, among which 42 are probable 
precursors of high-mass stars. This unbiased sample of young 
embedded massive stars reveals 17 dense cores that are ex- 
cellent candidates for harboring high-mass protostars in their 
infrared-quiet phase. The present paper also gives the first sta- 
tistical results on the lifetime of high-mass protostars and pre- 
stellar cores. Finally, Sect. |6] summarizes our conclusions, and 
Appendices A and B describe our data reduction and source ex- 
traction techniques. 

2. Observations and data reduction 

2.1. imaging the entire moiecuiar compiex 

2 

The Cygnus X molecular cloud complex covers up to ~ 30° 
(see Fig. 1 of Schneider et al. 2006) but the sites of high- 
mass star formation are found in its highest-density parts. 
We therefore used an infrared extinction map produced from 
2MASSC]data (cf. Fig. [T] and Bontemps et al. in prep.) to select 
the high-column density (Ay > 15 mag) clouds of Cygnus X 



' The Two Micron All Sky Survey was performed by the University 
of Massachussets and IPAC/Caltech and funded by NASA and NSF. 
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to be mapped in millimeter continuum emission. Figure [T] and 
[Schneider et al. (2006)| show that the Cygnus X molecular com- 
plex contains several ensembles of dense molecular clouds 
which extend over ~ 50 parsecs. In the following, we use the 
name "CygX-North" for the north-eastern group of clouds which 
harbors the well-known sources DR21 and W75N, "CygX- 
South" for the southern region containing AFGL 2591 and 
S106-IR, and "CygX-NW" for the north-western filamentary 
clouds. The mappings necessary to cover the high-column den- 
sity clouds of Cygnus X are outlined in Fig. [T] and cover a total 
area of ~ 3°'. 



observing runs was between 110 and 130 K and the average rms 
for each data set is 0.03 K on a main beam brightness temper- 
ature scale, taking an efficiency 77MB of 0.78. The angular reso- 
lution of the telescope at 87 GHz is 29". All observations were 
performed in position-switching mode with off-positions a few 
arcminutes away in azimuth. Pointing was checked regularly and 
was accurate to within 3". Standard calibration sources were ob- 
served regularly and were consistent within 10%. The data were 
reduced with the IRAM software for spectral lines (CLASS), ap- 
plying first order baselines and averaging individual spectra with 
rms weighting. 



2.2. Dust continuum observations at 1.2 mm 

The Cygnus X molecular complex was imaged at 1.2 mm 
with the MPIfR bolometer arrays installed at the IRAIvfl 30 m 
telescope at Pico Veleta (Spain). The passband of the MPIfR 
bolometers has an equivalent width of a; 70 GHz and is centered 
at Veff ~ 240 GHz (Krey sa et al.T99 8|. The MAMBO (37 chan- 
nels) camera was used in February-April 1999, January-March 
2000, and December 2000-May 2001 for a total integration time 
of ~ 30 hours. The MAMBO-2 (117 channels) camera was used 
for ~ 33 hours as part of the observing pools organized by IRAM 
between November 2002 and March 2003. A total of 142 large 
on-the-fly maps were taken using a fast-mapping technique sim- 
ilar to that described by Teyssier & Sievers (1998) These maps 
consist of a series of rows scanned in azimuth with a velocity of 
8"/sec and spaced by 22" (respectively 46") steps in elevation 
when using MAMBO or MAMBO-2. The typical azimuthal size 
of individual maps is ~ 7.5' with MAMBO versus ~ 11' with 
MAMBO-2. The signal of each bolometer is modulated by the 
secondary mirror which is wobbling with a frequency of 2 Hz 
and a throw in azimuth of 70". The resulting dual-beam maps 
were reduced with the IRAM software for bolometer-array data 
(NIC; cf.|Broguiere et al. 1995| l using the EKH restoration algo- 
rithm (lEmerson et al. 19791 1 and a skynoise reduction technique 
developed by us (see Appendix A). The imaging of each re- 
quested field was obtained by combining partially overlapping 
restored on-the-fly maps. 

The main beam was measured and found to have a HPBW 
size of ~ 11" using Uranus and Mars. The absolute pointing of 
the telescope was found to be accurate to within ~ 5". The data 
were all taken in winter but with heterogeneous weather condi- 
tions. The zenith atmospheric optical depth varied between 0.1 
and 0.5 and the skynoise level was low to high (correlation fac- 
tor from 0.2 to 1). The resulting rms noise in individual maps 
is cr = 10-150 mJybeam"', reduced to cr ~ 20 mJybeam"' 
after reduction of the skynoise (see Appendix A). The final mo- 
saics all have cr < 20 mJybeam"' rms with a median rms of 
cr ~ 15 mJybeam"'. Uranus and Mars were also used for flux 
calibration and the overall absolute calibration uncertainty is es- 
timated to be ~ 20%. 



2.3. SiO(2-1) observations 

Pointed observations in the SiO v=0 J=2— >1 transition were per- 
formed for 40 millimeter continuum sources listed in Table [T]in 
September 2003 and July 2004 with the IRAM 30 m telescope. 
We used the A 100 (BlOO) SIS receiver and the VESPA autocor- 
relator with a frequency resolution of 40 kHz that is equivalent 
to 0.135 kms"' at 86.847 GHz. The system temperature for both 



2 IRAM is supported by INSU/CNRS (France), MPG (Germany) and 
ION (Spain). 



3. Analysis of the 1.2 mm mapping of Cygnus X 

3. 1 . Dust continuum images 

The 1.2 mm continuum images are presented in Figs.|2^-c with 
zooms in the left-hand parts of Figs. [3]B1 (see also Figs. [T2HT41 
which are only available electronically). The dense molecular 
gas sometimes follows filamentary structures (well observed in 
the '^CO(l-O) images of Schneider et al. in prep.) and displays 
many compact (~ 0.1 pc) fragments. We make a complete cen- 
sus of these 1 .2 mm fragments in Sect. 13.21 and a similar but 
less complete census of larger-scale (~ 1 pc) cloud structures in 
Sect.|33] 

The MAMBO-2 mosaics are sensitive to spatial scales rang- 
ing from ~ 0.09 pc (corresponding to an HPBW of 11" at 
1.7 kpc) to ~ 5 pc (i.e. the mean azimuth extent of individual 
maps ~ 10' at 1.7 kpc). The resulting spatial dynamic is ~ 55, 
larger than that of other (sub)millimeter dust continuum obser- 
vations made with SCUBA at the JCMT, SIMBA at the SEST, 
or BOLOCAM at the CSO. Therefore, large fast-mapping im- 
ages with MAMBO-2 are currently the best tool to study the 
density structure of molecular clouds with high spatial resolu- 
tion. Given the dynamic range in the MAMBO-2 image (within 
all" beam, cold dust emission ranges from cr ~ 1 5 mJy beam" ' 
to 6 300 mJybeam"'), the Cygnus X complex is probed from 
A^H, ~ 9 X 10^' cm.-^ to 2 X lO^'* cm ^, or similarly from 
Ay ~ 10 mag to 2100 mag (i.e. with a column density dy- 
namic of more than 200). The above estimate uses Eq. 1 from 
[Motte et al. (I998)| with dust temperatures of 20 K and dust 
opacities of K12 mm - 0.005 and 0.01 cm^ g~' for the lower and 
higher density cloud structures, respectively. 

Due to the spatial filtering of cloud structures larger than 
clumps (i.e. > 5 pc, see Appendix A) in MAMBO-2 images, 
the total mass detected in millimeter continuum emission is 10 
to 20 times smaller (37 x 10^ Mq in CygX-North and 26 x 
10"^ Mq in CygX-South) than that measured in the CO survey 
of [Schneider etal. (2006)| 

3.2. Census of compact cloud fragments in Cygnus X 

Since our main goal is to investigate the best potential sites 
of high-mass star formation, we first focus on the small-scale 
cloud fragments observed in our MAMBO-2 images. As ex- 
pected, larger-scale (~ 1 pc) cloud structures generally have 
a moderate volume-averaged density (~ 5 x lO"' cm"-' for a 
1.2 mm peak flux of 5cr - 75 mJybeam"' over a diame- 
ter of 1 pc, when a 20 K dust temperature and a /ci jmm = 
0.01 cm^g"' dust mass opacity are assumed, cf. Eqs. [T] and 
O and will therefore be less likely to form high-mass stars in 
the near future. We apply the source extraction technique devel- 
oped by Motte et al. (2003) which uses a multi-resolution anal- 
ysis ( [Starck & Murtagh 2006| ) and the Gaussclumps program 
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Fig. 2. Millimeter continuum imaging of the Cygnus X molecular cloud complex obtained with the MAMBO and MAMBO-2 
cameras installed at the IRAM 30 m telescope. These 1.2 mm maps have been smoothed to an effective angular resolution of 15", 
allowing a sensitivity of 0. 1-5 pc cloud structures. The main radio sources (Downe s & Rinehart 1 966 ) and a few well-known sources 
are indicated as reference marks, a The CygX-North region (see Fig.[T]for its location): maximum flux is ~ 8 500 mJy beam"' (color 
scale is saturated beyond 500 mJy beam"') and rms noise level is cr = 10-20 mJybeam"'. 
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Fig. 2. (continued) b The CygX-South region: 



maximum flux 

-1 



is ~ 8 500 mJybeam (color scale is saturated beyond 



500 mJybeam" ) and rms noise levels are 9 - 18 mJybeam" in bl, ~ 14 mJybeam" in b2, ~ 14 mJybeam" in b3, 
~ 13 mJybeam ' in b4, and ~ 15 mJybeam ' in b5. Note that the MAMBO-2 maps associated with DR12 and DR6 in Fig.[T] 
are ignored since they do not display any compact and dense gas. 
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(c) CygX-NW region at 1.2 mm 




Fig. 2. (continued) c The CygX-NW region: maximum flux 
is ~ 800 mJybeam"' (color scale is saturated beyond 
500 mJybeam"') and rms noise levels are 9-16 mJybeam"' 
in cl, ~ 13 mJybeam ' in c2, ~ 18 mJybeam ' in c3, and 
~ 12 mJybeam"' in c4. 



dStutzki & Giisten 19901 IKrameretal. 1998l l. We first identify 
129 compact (< 1 pc) cloud fragments above the median 5cr - 
75 mJybeam"' level in the MAMBO-2 maps, where spatial 
scales larger than 1 pc have been filtered out (see e.g. Fig.fTTb). 
We then precisely derive the 1 .2 mm characteristics of each frag- 
ment by fitting a 2D-Gaussian in the image where its local back- 
ground (or surrounding large-scale emission) has been entirely 
removed (see Appendix B for more details). 

The compact fragments extracted here are marked and la- 
beled in the left-hand parts of Figs. [3}E] and Figs. \T2WT4\ 
Table [T] which is only available electronically, lists the 129 
fragments discovered in Cygnus X and gives their 1 .2 mm ba- 
sic characteristics as estimated above. The fragments are listed 
in increasing RA coordinate for each region. They are num- 
bered (Col. 1) and named after their 2000 equatorial coordi- 
nates (Col. 2). Their 1.2 mm peak flux density (Col. 3), their 
FWHM size (Col. 4), and integrated flux (Col. 5) are mea- 
sured assuming a Gaussian shape. In Col. 9, we list previous 
names given to strong millimeter peaks of the DR21 filament 
(cf. lChandler et al. 1993al l. as well as the name of nearby bench- 
mark sources (e.g. [Carey et al. 2000[ ). 



3.3. Large-scale cloud structures in Cygnus X 

To complement the study performed in Sect. l3.2l we here make a 
census of the large-scale (~ 1 pc) cloud structures of Cygnus X. 
We first smooth the MAMBO-2 maps to a 55" resolution, which 
roughly corresponds to ~ 0.5 pc and is the resolution of a N2H^ 
survey performed at the FCRAO (Schneider et al. in prep.). We 
then identify 40 large-scale cloud fragments above the median 
5cr = 75 mJy beam"' level and fit them by a 2-D Gaussian. Note 
that the skynoise reduction technique we applied may have fil- 
tered out part of the large-scale emission (for more details see 
Appendix A). The present analysis is therefore less secure than 
that performed in Sect. l3.2l for the small-scale fragments of the 
Cygnus X molecular clouds. Table |2] which is only available 
electronically, lists the 40 large-scale clumps found in Cygnus X 
and gives their 1.2 mm basic characteristics (in Cols. 1-6) with 
the same convention as in Table [T] In addition. Col. 7 lists the 
compact fragments (if any) harbored by each larger-scale struc- 
ture. 

3.4. Mass and density derived from millimeter continuum 
maps 

The mass and density of the compact cloud fragments (< 1 pc) 
and large-scale structures (~ 1 pc) of Cygnus X are determined 
from their dust continuum parameters given in Tables [T] and |2] 
These tables list the mass (Col. 6 vs 5) and volume-averaged 
density (Col. 7 vs 6) estimated for each identified dense core 
and clump using methods and assumptions described below. 

Before we estimate the mass of each cloud fragment, 
we should ensure that the millimeter sources identifed in 
Sects. 13.21 and 13.31 are indeed density structures. According 
to previous studies of high-mass star forming regions (see 
e.g. Mot te et al. 20031 1, the main source of contamination of 
1.2 mm fluxes is the free-free emission from H ii regions. 
Unfortunately, the published centimeter wavelength free-free 
images of Cygnus X (e.g. IWendker et al. 199l1 l lack the neces- 
sary angular resolution to estimate the amount of likely contami- 
nation of each 1 .2 mm source. Furthermore, the spectral index of 
the free-free emission is unknown, precluding the extrapolation 
of its flux at 1 .2 mm. Observations with high-density tracers such 
as molecular lines at the IRAM 30 m, Eff'elsberg, and FCRAO 
telescopes, as well as 350 yum continuum with SHARC II at the 
CSO, suggest that the vast majority of the 1.2 mm compact and 
larger-scale fragments are definitive density structures traced 
by dust emission rather than free-free centimeter sources (e.g. 
Schneider et al. in prep.). Fifteen millimeter sources of Table [T] 
are suspected to contain H n regions (cf. Sect. l5.2.TT i; their asso- 
ciated masses and densities are therefore upper limits. 

We believe that the millimeter continuum emission of cloud 
fragments identified in Sects. [l!2] and l3.3l is mainly thermal dust 
emission, which is largely optically thin. For any given dust 
properties and gas-to-dust ratio, the 1.2 mm fluxes are thus di- 
rectly related to the total (gas -i- dust) mass of the fragments. 
For present mass estimates, we use the 1 .2 mm integrated fluxes 
(S j'"'j^j^) of Tables[TJ[2]without correction for any free-free con- 
tamination, and derive the mass (Mi 2 mm^ given in Tables [T]-|2]i, 
as follows: 

5/"' 

_ 1-2 mm 

^"1.2 mm — n (T \ 

1^1.2 mm ^1.2 mml-' dust/ 
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1^1.2 mm 

0.01 cm^g- 



\ ' /£dust\ 
"/ \20K/ 



(1) 



The volume-averaged densities listed in Tables[T]-(2]and used 
in Tables |4]-[5] are estimated as follows: 



where /cj 2 mm is the dust opacity per unit mass column density 
at 1.2 mm, B[ 2 mm(7^dust) is the Planck function for a dust tem- 
perature Tjust, and where we assume a distance of 1.7 kpc. 

The dust mass opacity (including dust properties and gas- 
to-dust mass ratio) is likely to vary with density, temper- 
ature, and the evolutionary state of the emitting medium 
( [Henning et al. 1995). Models of dust i n low-mass protostellar 
cores (e.g. |Ossenkopf & Henning 1994 1 suggest that a value of 
''■1.2 mm - 0.01 cm^g"' is well suited for cool (10 - 30 K) 
and high-density (hht ~ 10^ cm"^) cloud fragments. This value 
agrees with the cross-comparisons of dust emission surrounding 
an UCH 11 region with its CO ice absorption and gas emission 
dvan der Tak 2002t . We therefore choose a dust opacity per unit 
(gas + dust) mass column density of '<'i.2mm = 0.01 cm^g"' 
for compact cloud fragments. We use this value for all the dense 
fragments because the average dust properties of these ~ 0. 1 pc 
structures should not change drastically when they contain a 
massive protostar. We also adopt /ci.2mm = 0.01 cm^g"' for 
larger-scale (~ 1 pc) clumps, despite the fact that a lower value 
is generally used for lower-density starless cloud fragments (e.g. 
Ward-Thompson et al. 1999). We estimate that the absolute value 
taken for the dust mass opacity is uncertain by a factor of 2. 

The temperature to be used in Eq. ([T) is the mass-weighted 
dust temperature of the cloud fragments, whose value could be 
determined from gray-body fitting of their spectral energy dis- 
tributions. Such measures are not yet available but we recently 
measured the ammonia rotational temperatures of the Cygnus X 
dense cores (Wyrowski et al. in prep.). The ratios of NH3 (1,1) 
and (2,2) transitions observed with the Effelsberg 100 m tele- 
scope (40" beam) give cold temperatures (~ 15 K) for the ma- 
terial located within a 0.33 pc diameter of a millimeter peak. 
To estimate the temperature averaged on ~ 0.1 pc scales, we 
assume a classical p(r) oc r^^ density law and inner heating 
(T{r) oc r""-^), in agreement with small-scale fragments be- 
ing highly centrally condensed and generally protostellar (cf. 
Sects. I5TI45.2I 1. We thus estimate that ~ 0.1 pc fragments should 
have mass- weighted dust temperatures in the range of 15-25 K. 
This does not preclude any hotter or colder components from 
being embedded in the ~ 0.1 pc fragments; temperature gra- 
dients will be discussed in companion papers (Schneider et al. 
in prep.; Wyrowski et al. in prep.). For simplicity, we assume 
^dust - 20 K in Eq. ([TJ for all the cloud fragments of Table [1] 
We add estimates made with Tjust = 40 K for H n regions 
(identified in Sect. 15.2. il l and the brightest (Lbol > 5 x IO^Lq) 
infrared protostars of C ygnus X: IRAS 20343+4129, studied 
by Beut her et al. (2002)] which corresponds to CygX-N5, and 
LkHA 225S which coincides with CygX-NW5. To avoid bias- 
ing our statistical results, we hereafter use the estimates made 
with a homogeneous 20 K temperature. The dust temperature of 
large-scale clumps is more difficult to constrain since it may vary 
according to the presence or absence of strong external heating. 
We estimate that the mass-weighted dust temperature may vary 
from r^just = 10 K to 20 K, and choose Tjust = 15 K as an aver- 
aged value to use in Eq. ([TJ for the large-scale (~ 1 pc) clumps 
of Table 121 

The mass estimates in Tables [T]-(2] are globally uncertain by 
a factor of 2 due to uncertain dust emissivity. The individual 
masses can vary by ±30% relative to each other when dust tem- 
peratures vary from 15 to 25 K, or by +50% for the 10 to 20 K 
temperature range (see also Eq.[T]). 



; FWHM' ' 



(2) 



where Mi 2mm is the mass derived by Eq. ([T) and FWHM is 
the full width at half maximum, determined by Gaussian fits. 
Using a radius equal to one FWHM in Eq. (|2]i allows one to 
accurately determine the volume-averaged density because the 
flux (and thus the mass) measured within such an aperture corre- 
sponds to > 98% of the integrated flux (respectively total mass) 
of Gaussian cloud structures. The often-used, beam-averaged 
peak density would of course be higher, but is less relevant when 
estimating physical constraints such as the free-fall time. 

4. Signpost of stellar activity toward the Cygnus X 
cloud fragments 

We investigate the spatial coincidence of the compact cloud frag- 
ments identified in Sect. |3.2| with mid-infrared, SiO, centimeter 
free-free, and maser emission (see Sects. 14.1144.31 1. Large-scale 
cloud structures identified in Sect. 13. 31 which are associated with 
neither compact fragments displaying signposts of stellar activ- 
ity nor bright mid-infrared sources, qualify as starless clumps as 
indicated in Col. 7 of Table |2] 

4.1. Coincidence witli mid-infrared point sources 

When searching for the signature of embedded stellar embryos, 
the cross-correlation of compact millimeter sources with in- 
frared catalogs is crucial. Among the near- to mid-infrared sur- 
veys that are currently available and complete for the Cygnus X 
complex, we choose that performed by the satellite; its 

angular resolution (20") is far better than that provided by 
IRAS and its wavebands (from 8 to 21 fim) are better suited 
than those of 2MASS. The Spitzer satellite provides better an- 
gular resolution and sensitivity but the mapping of Cygnus X 
with the Spitzer/MlPS camera is far from being complete, 
and published results from the Spitzer/IRAC camera generally 
do not cross-correlate well with (sub)millimeter sources (see 
Davis et al. 2007 and the counterexample ER01/CygX-N57 of 
Marston et al. 200 4). 

To search for coincidence with mid-infrared point sources, 
we carefully examine the structure of the MSX emission at the 
exact location of each compact millimeter source (see right-hand 
side of Figs. [3]-@]and Figs. [l"2l - [T4l i. A millimeter source is de- 
fined to be detected by MSX if it coincides with pointlike emis- 
sion at both 8 and 21 fim. Slightly more extended MSX emission 
is considered for three developed H 11 regions (cf. Sect. 15.2. Il l: 
DR21 (associated with CygX-N46 and N47), S106-IR (associ- 
ated with CygX-S18, S19 and S20) and IRAS 20306+4005 (as- 
sociated with CygX-S39 and S40). Given the resolution of MSX 
images, MSX point sources have similar sizes to those of com- 
pact cloud fragments detected by MAMBO-2. The detection, at 
both 8 and 21 yum, allows the rejection of peaks of infrared emis- 
sion created by small grains at the interface of photo-dissociation 
regions and molecular clumps (e.g. Abergel et al. 2002| ). The 
latter are actually externally heated sources whose peaks shift 
with wavelength and sometimes disappear at 21 pm. We add 
four more sources which are only detected at 21 pm: CygX- 
N30 (W75N(B)), CygX-N46 (DR21) and N47 (DR21-D), whose 



^ The Midcourse Space Experiment provides images observed at 8, 
12, 15, and 21 //m with a 20" resolution ^Egan et al. 2001| l. 
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Fig. 3. MAMBO maps of CygX-North (left: gray-scale and contours, right: contours overlaid on 8 yum images obtained by MSX 
and converted to Jy sr ') extracted from Fig. |2^. Regions shown are west of DR21 (a), around DR21 (b) and DR22 (c). The 
1.2 mm and 8 images have 11" and 20" angular resolutions, respectively. The compact cloud fragments discovered in MAMBO 
images (see Table[T]i are labeled and marked by crosses in the gray-scale plot. The infrared sources which coincide with a MAMBO 
cloud fragment are also indicated. Contour levels are logarithmic and go from 40 to 800 mJy beam"' in a and c, and from 40 to 
4 800mJybeam ' in b. 
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Fig. 4. Same convention as Fig.[3]for MAMBO maps of CygX-South extracted from Fig.|2]3. Regions shown are south of DR4 (a), 
around S106 (b), south of DR15 (c), and around DR15 (d). Contour levels are logarithmic and go from 40 to 200 mJy beam ' in a 
and d, from 60 to 1 200 mJy beam ' in b, and from 40 to 400 mJy beam ' in c. 
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8 //m emission is extinguished by their high-column density, 
and CygX-N65, which is very weak. The positional accuracy re- 
quested between the compact sources of MAMBO-2 and MSX 
is 10" (i.e. the sum of the maximum pointing errors of both 
images), which corresponds to 0.08 pc at 1.7 kpc. Such a nar- 
row association is made possible by the small size (~ 0.1 pc) 
of the compact cloud fragments extracted from the MAMBO-2 
images. In fact if we consider looser associations, many MSX 
point sources lie on and possibly within clumps, at ~ 0.2 pc 
from the density peaks which we identified as compact frag- 
ments (e.g. Figs. [3^-b). We reject these loose associations, un- 
like previous studies that considered, for example, that CygX- 
S12 (respectively CygX-S30) and IRAS 20205+3948 (respec- 
tively IRAS 20293-H3952) coincide (cf. Figs. |4}i and Ej;, see 
also Beuther et al. 2002). In fact, these MSX point sources ly- 
ing within dust clumps sometimes coincide with extremely red 
objects detected by Spitzer/IRAC (EROs, iMarston et al. 20041 
IDavis et al. 20071) . confirming that (sub)milUmeter and infrared 
surveys generally identify different population of stars in terms 
of age and/or mass. 

Our cross-correlations are illustrated in the right-hand parts 
of Figs. [3]B1 and Figs. [T2HT4l Their results are given in Col. 8 of 
Table [1] When there is a coincidence with a MSX point source. 
Col. 8 gives its 21 fim flux as taken from the MSX catalog (see 
Note c for exceptions) and Col. 9 gives the IRAS or MSX name. 
Otherwise, when the MAMBO-2 source is not detected or seen 
in absorption at 8 /vm. Col. 8 of Table [T]mentions "-" or "Abs". 

With the above selection criteria, 36 compact cloud frag- 
ments of Cygnus X (i.e. less than 1/3 of the sample) are as- 
sociated with mid-infrared emission detected by MSX. The other 
93 fragments are undetected down to ~ 0.15 Jy at 8 fim and 
~ 3 Jy at 21 //m, which are the completeness levels of the 
MSX catalog (MSX C6) in Cygnus X (see Egan et al. 200T| 



Bontemps et al. in prep.) . Some of these sources are even seen 



in absorption against the diff'use mid-infrared background. 
Interestingly, some of the Cygnus X fragments which are un- 
detected by MSX harbor a weak 24 fim source in Spitzer images. 
To preserve the homogeneity of our study, we do not use this 
information here. 



4.2. Survey for SiO emission 

Silicon monoxide (SiO) emission is an excellent tracer 
of shocked gas, usually associated with molecular out- 
Ziurys & Friberg 1987[ IMartm-Pintado et al. 19921 



flows (e.g. 

ISchilke et al. 1997|). SiO may also have been detected in 



shocks toward hot cores (e.g. IHatchell et al. 200T 1 and refer- 
ences therein). We have surveyed for SiO(2 - 1) emission by 
making single pointings toward almosQ all the most massive 
(Ml 2 mm > 40 Mq) fragments of Table[T]and a few less massive 
fragments. Figure |5] and Fig. [15] (only provided electronically) 
display the SiO(2 - 1) lines for all the observed fragments, or- 
dered by decreasing Mj 2mm- In Tableland for each 1.2 mm 
compact fragment (Col. 1), we give the peak line temperature 
(Col. 3), the line width at half maximum (Col. 4), the full line 
width at the base (Col. 5), and the integrated intensity (Col. 6) 
of the SiO(2-l) emission. For direct comparison, we list the SiO 
line parameters scaled to the Cygnus X distance of three nearby 



* Three sources (CygX-N32, N47, and S20) are included in beams 
targeting other fragments (CygX-N30, N46, and S 1 8 respectively). Two 
sources (CygX-N69 and CygX-SlO) were not observed because they 
are among the lowest-density fragments. 



protostars which are among the brightest known SiO outflow 
sources. 

Among the 28 most massive fragments whose SiO emis- 
sion is displayed in Fig. |5] only two (CygX-S18, and CygX- 
N6) are not detected. This implies a very high detection rate of 
93 % (> 79 % when considering the five high-mass fragments 
which have not been observed), with a typical detection level 
of ~ 0.15 Kkms ' (see Table|3ll. This detection rate stays high 
even when considering the observation of the 20 - 40 Mq dense 
cores listed in Table [3] Down to ~ 25 Mq, the compact frag- 
ments of Table[3]that do not have any SiO emission are all asso- 
ciated with infrared point sources: CygX-S18/S106-IR, CygX- 
N6/IRAS 20343+4129, and CygX-N16/MSX 81.4452+0.7635. 
This suggests that all massive fragments down to at least 40 Mq, 
and maybe 25 Mq, drive SiO outflows and/or contain hot cores. 

In Fig. |5] we show the emission spectra separately for the 
high-luminosity infrared sources (> 10^ Lq) and the massive 
infrared-quiet cores (cf. their definitions given in Sects. 15^2721 and 
15.3.1b . The two groups obviously have different behavior. All 
infrared-quiet cores are detected in SiO(2-l) and most of them 
are particulary bright; the averaged peak line temperature and 
integrated intensity are 0.55 K and 4.6 Kkm s"', respectively. In 
contrast, only a small fraction of the high-luminosity infrared 
sources show a strong SiO line; the averaged peak line tem- 
perature and integrated intensity are 0.22 K and 1.9 Kkms"', 
respectively. Since only two high-luminosity infrared sources 
(CygX-N48/DR21 and CygX-N30/W75N) are strong SiO emit- 
ters, the contrast is even larger when comparing the median 
values: 0.7 Kkms"' for the high-luminosity infrared sources, 
and 3.1 Kkms"' (i.e. 4.3 times more) for the infrared-quiet 
cores. The SiO lines of infrared-quiet cores are also 3 to 4 times 
brighter than the brightest SiO outflows of the nearby low-mass 
protostars (see Table[3]l. 

Almost all the SiO lines shown in Figs. |5] and [15] dis- 
play outflow wings with line width at the base of emission 
that can reach up to ~ 60 kms"'. Among the ten brightest 
SiO sources, the average of the line width at the base is as 
high as 41 kms"', thus clearly suggesting outflow shocks. A 
few SiO lines might also indicate the presence of at least two 
emission components: a wide base (due to outflow shocks) 
and a narrower component which may indicate a contribu- 
tion of shocks inside a hot core and/or lower velocity out- 
flows. The best examples of two-component lines are those 
detected toward CygX-S26/AFGL 259, CygX-N44/DR21(OH), 
and CygX-N48/DR21(OH)-S, while the best example of single- 
component lines suggestive of a unique powerful outflow are 
CygX-N46/DR21, CygX-N53, CygX-S30, and CygX-N68. 

4.3. Coincidence witli radio free-free and maser emission 

We have used the SIMBAEjfl database to search for complemen- 
tary signposts of stellar activity: centimeter free-free emission 
and OH, H2O, and CH3OH masers. SIMBAD gives inhomoge- 
neous information on our sample of 129 compact cloud frag- 
ments, since it compiles results of Cygnus X surveys performed 
with limited angular resolution and high-angular resolution stud- 
ies made for well-known sources. 

Within Table [T] fifteen dense cores coincide with a 
small-diameter (< 10") source of radio emission de- 
tected at 2, 3.6, 6, 11, or 21 cm (e.g. IKurtz et al. T994t 



' SIMBAD is the reference database for identification and bibli- 
ography of astronomical objects, developed and maintained by CDS, 
Strasbourg. 
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Fig. 5. SiO(2-l) lines observed toward the 28 most massive (> 40 Mq) cloud fragments of Cygnus X: 14 high-luminosity infrared 
sources on the left and 14 infrared-quiet protostellar cores on the right, both ordered by decreasing My 2 mm from top to bottom. 
The local velocity at rest of each cloud fragment, as measured by optically-thin tracers, is indicated by a dashed line. 



IZoonematkermani et al. 199 ^. We cannot check the na- 
ture of the detected centimeter emissions with the cur- 
rent database, but it most likely corresponds to free- 
free emissions emitted by H 11 regions. Furthermore, sev- 
eral 1.2 mm continuum sources coincide with OH, H2O, 
and/or CH3OH mase r sources (e.g. |Braz & Epchtein 1983[ 
IValdettaro et al. 20011 IPestalozzi et al. 20051) . Most of these 



dense cores are associated with centimeter continuum and/or 
mid-infrared emission, but one (CygX-N53) is neither a cen- 
timeter nor a MSX source. In Col. 9 of Table[T] we indicate when 
some centimeter free-free and/or maser emission is detected to- 
ward a Cygnus X dense core. 
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Table 3. SiO(2-l) observations of the most massive dense cores of Cygnus X compared with the brightest SiO protostellar sources 
at < 500 pc 



Source 


A?1.2mm 


7'sio 


A VfWHM 


^ ^'outflow 


J Tsio dv 


Possible nature 


name 


(Mo) 


(K) 


(km s"') 


(km s"') 


(Kkms"') 




N46, DR21 


< 609 - 949 


0.38 


16.9 


64.8 


8.6 ± 0.42 


UCH II region 


N30, W75N(B) 


< 361 - 563 


0.65 


5.4 


25.8 


4.7 ±0.16 


hyper compact H ii region? 


N44, DR21(OH) 


446 


1.45 


6.0 


43.0 


13.4 + 0.25 


massive IR-quiet protostellar core 


S18, S106-1R 


< 145 - 226 


< 0.10 


- 


- 


<0.15 


compact H ii region 


N48, DR21(OH)-S 


197 


0.78 


4.0 


49.7 


5.4 ± 0.26 


massive IR-quiet protostellar core 


S8, Mol 121S 


< 121 - 188 


0.10 


4.6 


22.1 


0.96 + 0.18 


UCH II region 


N58, IRAS 20375+4109 


< 92 - 144 


0.086 


3.8 


7.1 


0.34 + 0.10 


compact H ii region 


N38, DR21(OH)-W 


138 


0.86 


4.6 


37.5 


6.8 + 0.25 


massive IR-quiet protostellar core 


N51, W75S-F1R2 


125 


0.29 


3.6 


12.7 


1.5 + 0.19 


high-luminosity IR protostellar core 


S26, AFGL2591 


< 80 - 124 


0.30 


2.8 


18.7 


1.5 + 0.19 


UCH II region 


N40, DR21(OH)-N2 


106 


0.63 


2.6 


14.4 


2.3 + 0.18 


massive IR-quiet protostellar core 


N12 


86 


0.87 


4.0 


26.5 


4.7 + 0.21 


massive IR-quiet protostellar core 


N53 


85 


0.80 


6.4 


57.6 


8.9 + 0.31 


massive IR-quiet protostellar core 


N3 


84 


0.40 


4.8 


21.5 


2.6 + 0.21 


massive IR-quiet protostellar core 


S30 


82 


0.32 


9.3 


35.2 


3.9 ± 0.26 


massive IR-quiet protostellar core 


NIO, IRAS 20350+4126 


<52-81 


0.12 


2.4 


4.7 


0.35+0.10 


UCH II region? 


N6, IRAS 20343+4129 


44-68 


< 0.07 


- 


- 


<0.10 


high-luminosity IR protostellar core 


N43, W75S-FIR1 


66 


0.25 


3.4 


9.1 


0.94 + 0.15 


high-luminosity IR protostellar core 


NWl,IRAS20178+4046 


< 38-59 


0.05 


2.6 


3.4 


0.15+0.06 


UCH II region? 


N63 


58 


0.23 


7.2 


36.5 


2.6 ± 0.23 


massive IR-quiet protostellar core 


S32 


54 


0.29 


6.2 


24.6 


2.4 ± 0.28 


massive IR-quiet protostellar core 


N14, IRAS 20352+4124 


50 


0.16 


2.4 


7.1 


0.61 +0.13 


high-luminosity IR protostellar core 


NW14, Mol 124 


32-50 


0.16 


41.6 


64.4 


5.8 ± 0.76 


high-luminosity IR protostellar core 


S37, IRAS 20305+4010 


45 


0.05 


- 


3.2 


0.09 ± 0.09 


massive IR-quiet protostellar core 


N68 


44 


0.27 


13.7 


43.8 


4.6 ± 0.25 


massive IR-quiet protostellar core 


S41, IRAS 20306+4005 


27-42 


0.30 


2.0 


8.2 


0.61 +0.09 


high-luminosity IR protostellar core 


N65 


41 


0.19 


3.8 


16.9 


1.06 + 0.16 


massive IR-quiet protostellar core 


S7 


40 


0.58 


6.0 


28.4 


5.1+0.23 


massive IR-quiet protostellar core 


NW2 


36 


0.10 


- 


2.0 


0.18 + 0.09 


IR-quiet protostellar core 


S43 


35 


0.16 


2.2 


6.1 


0.59 + 0.13 


IR-quiet protostellar core 


N24 


34 


0.09 


4.7 


5.1 


0.39 + 0.10 


IR-quiet protostellar core 


NW5, LkHA 225S 


22-34 


0.045 


- 


14.9 


0.69 ± 0.26 


infrared protostellar core 


N60 


30 


0.08 


2.0 


4.5 


0.20 ± 0.10 


IR-quiet protostellar core 


NI6, MSX 81.4452+0.7635 


29 


< 0.09 






< 0.12 


infrared protostellar core 


INzo 


L 1 


u.uy 


1 o 
1.0 


J.J 


U.Zj ± U.uy 


IR-quiet protostellar core 


N62 


27 


0.10 


2.6 


8.5 


0.42 + 0.11 


IR-quiet protostellar core 


SI6 


24 


0.10 


2.6 


14.5 


0.51+0.13 


IR-quiet protostellar core 


NI8 


23 


< 0.08 






< 0.11 


? 


N59 


20 


< 0.08 






< 0.11 


? 


NW9, IRAS 20216+4107 


19 


< 0.10 






<0.15 


infrared protostellar core 


L1448-MM @ 1.7 kpc 








~ 140 


~ 1.4*"' 


low-mass class protostar 


L1157-MM @ 1.7 kpc 








-25 


~ 1.8<''> 


low-mass class protostar 


Orion-IRc2 @ 1.7 kpc 








~ 74 


~ 9.2<=' 


high-mass protostar 



" Intensity which has been distance-corrected (extrapolation from 300 pc to 1.7 kpc) using values from Table 2 of |Nisini et al. (2007)| and 
excluding R4 which would fall outside the beam. 

Intensity which has been distance-corrected (extrapolation from 440 pc to 1.7 kpc) using values from Table 3 of [Nisini et al. (2007)| and 
summing up BO and B 1 components. 

Intensity which has been distance-corrected (extrapolation from 450 pc to 1.7 kpc) using the spectrum of Fig. 1 in Ziurys & Friberg (1987). 



5. Discussion 

5.1. A sample of high-density cores 

Our 1.2 mm continuum imaging of the Cygnus X molecu- 
lar cloud complex gives a complete view of the structure of 
its highest-column density portions, with spatial scales ranging 
from 0.03 pc (deconvolved from the beam) to 3 pc. According 
to the terminology used for nearby molecular clouds (e.g. 
IWilliams e t al. 2000), it covers scales typical for dense cores 
(~ 0.1 pc) and clumps (~ 1 pc). The 129 compact cloud frag- 
ments identified in Cygnus X have sizes ranging from 0.03 to 



0.3 pc with a mean size of ~ 0.1 pc (see Sect. l3.2l and Table[Tll. 
They should therefore be called "dense cores". The 40 larger- 
scale structures of Table |2] have typical sizes of ~ 0.7 pc; 
we will call them "clumps". In Table |4] we give the mean 
properties of the Cygnus X dense cores and clumps and com- 
pare them with those of a few representative studies (Motte 
et al. 1998; Ward-Tho mpson et al. 1999; Beuther et al. 2002; 
Rat hborne et al. 2006l l. In Fig. |6^-b, we plot the mass and den- 
sity distributions of these 1.2 mm sources samples. 

The typical size of the Cygnus X dense cores is 5 to 10 times 
smaller than the average cloud structures identified by single- 
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Table 4. Mean properties of the Cygnus X clumps and cores compared with millimeter cloud structures of a few reference studies 





HMPOs 


IRDCs 


r^vQ^ pliiTTins 




> 40 CvpX 




p Oph 




clumps 


clumps 


(all / starless) 


dense cores 


dense cores 


dense cores 


condensations 


FWHM " size (pc) 


0.5 


0.5 


0.68/0.78 


0.10 


0.13 


0.08 


0.007 


Mass* Ml 2 mm(A^o) 


290 


150 


1000 / 780 


24 


91 


4.7 


0.15 




8.5 X 10' 


5.9 X 10' 


1.4x 10V6.7X 10' 


1.1 X 10^ 


1.9 X 10' 


3.5 X 10" 


1.9 X 10* 


Number of sources 


69'- 


190 


40/ 10 


128 


33 


22 


60 


References'' 


(1) 


(2) 


(3) 


(3) 


(3) 


(4), (5) 


(5) 


Distance (kpc) 


0.3-14 


1.8-7.1 


1.7 


1.7 


1.7 


0.14-0.44 


0.14 



(a) Deconvolved FWHM sizes derived from a 2D-Gaussian fit, except for HMPOs and IRDCs for which sizes are taken from Table 1 of 
[Beuther et a l. (2002) (mean estimates are taken when the distance is unknown) and Table 3 of Rathbome et al. (2006). 

(b) Mass consistently estimated from the 1.2 mm integrated flux using Eq. (1). The mass values of |Beuther et al. (2002)| have been recalculated 
using Ki 2 mm = 0.01 cm" g"' but those of nearby, low-mass dense cores and p Oph condensations are kept while they assume /cj 7 mm = 
0.005 cm^g-'. 

(c) Volume-averaged density homogeneously recalculated using Eq. ^ for every list of sources. 

(d) References: (l) |Beuther et al. (2002) 1 (2) Rathbome et al. (2006); (3) this paper; (4) [Ward-Thompson et al. (1999)[ (5) |Motte et al. (1998)| 

(e) We limit our comparison to the main HMPOs components, which are the most likely progenitors of high-mass stars. 




Fig. 6. Distribution of mass (a) and volume-averaged density (b) of the Cygnus X dense cores compared to those of clumps 
hosting high-luminosity protostars (main component of HMPOs, Beuther et al. 2002) and clumps referred to as IRDCs 
jRathborne et al. 2006 ). The Cygnus X histograms are also compared to density spectra of nearby dense cores (Motte et al. 1998; 
Ward-Thompson et al. 1999) and p Oph condensations (Motte et al. 1998). 



dish studies of high-mass star-forming regions (see Table H] and 
also Reid & Wilson 2005; Hill et al. 2005). Using the terminol- 
ogy mentioned above, HMPOs (for high-mass protostellar ob- 
jects, Beuther et al. 2002) and IRDCs (for infrared dark clouds, 
e.g. IRathborne et al. 20061 ) are actually "clumps". We are here 
able to reach the scale of dense cores because Cygnus X is one 
of the closest high-mass star-forming regions, our MAMBO-2 
survey has a good spatial resolution, and we deliberately fo- 
cus on high-density cloud structures. As a consequence of their 
small sizes, the Cygnus X dense cores have masses (4 - 950 M©, 
see Table [1) which are up to one order of magnitude smaller 
than those obtained by Beuther et al. (2002)| and Rathborne et 
al. (2006) (see Table |4j and Fig. |6Ji). In contrast, Cy gnus X 
dense cores have densities averaged over their FWHM size 
(1 X 10"* - 2 X 10^ cm"-') which are more than one order of 
magnitude higher than those of HMPOs and IRDCs (cf. Table |4] 
and Fig.jSJ?). As interesting illustrations, the Cygnus X HMPOs 
called IRAS 20343 -h4 129 and IRAS 20293 -h3952 (Beuther et 
al. 2002) are twice as large and 7 times less dense than the dense 
cores we have identified within them (see CygX-N5, CygX-N6 
and CygX-S30inTable[I]i. 



On average, the Cygnus X dense cores are a factor of 5 more 
massive and a factor of 3 denser than nearby dense cores iden- 
tified by Mo tte etal. (1998)| and [Ward-Thompson et al. (19"99) 
(see Tableland Figs.[6^-b). The latter sample includes isolated 
pre-stellar cores located at 140-440 pc from the Sun, and all the 
protocluster dense cores of the p Ophiuchi main cloud which are 
either pre-stellar or already forming low-mass stars. Statistically, 
studying the Cygnus X region gives access to dense cores which 
are one to two orders of magnitude more massive than those ob- 
served in nearby regions. The 33 most massive compact cloud 
fragments of Cygnus X (i.e. with masses larger than 40 Mg, see 
Table O already have 19 times more mass than the dense cores 
of low-mass star-forming regions (cf. Table |4|l. Thus, they do 
not have any equivalent in the nearby star-forming regions and 
represent good candidate sites for forming intermediate- to high- 
mass stars. 

In terms of average density, Cygnus X dense cores are in- 
termediary cloud structures between nearby dense cores and 
condensations in the p Ophiuchi main cloud (cf. Table |4| and 
Fig. I6J3). This result suggests that the efficiency of mass trans- 
fer from the gas reservoir to the star(s) is high. Indeed, the star 
formation efficiency estimated in the p Ophiuchi dense cores is 
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~ 31% over 10^ yr ( [Bontemps et al. 2001[ i, and the good agree- 
ment in p Oph between the pre-stellar condensation mass spec- 
trum and the stellar IMF suggests an efficiency of 50% - 100% 
(Motte et al. 1998; [Bontemps et al. 2001| l. While the low-mass 
starless condensations are believed to be the direct progenitors 
of single stars, the Cygnus X dense cores will probably form 
small groups of stars. The large mass of the Cygnus X dense 
cores, compared to that of nearby dense cores and p Oph conden- 
sations, argues for the formation of intermediate- to high-mass 
stars (cf. Fig. |6^). With their high density and small size, the 
Cygnus X dense cores can probably be seen as the inner part of 
clumps/protoclusters over which the stellar IMF does not glob- 
ally apply. Therefore, we may expect that the > 40 Mq dense 
cores of Cygnus X have a high probability of forming 10-20 
of stars, including at least one high-mass star. The detection 
of maser emission toward a few moderate-mass (~ 20 Mq) 
dense cores that remain weak at mid-infrared wavelength (cf. 
Wyrowski et al. in prep.) suggests that the > 40 M© limit for a 
Cygnus X dense core to host a high-mass protostar is reasonable. 



5.2. The high potential of Cygnus X dense cores to form 
massive stars 
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Fig. 7. Separating the high-luminosity (> 10 Lq) sources from 
infrared-quiet cores using their 21 fim flux (limit set to 10 Jy). In 
Cygnus X, 15 UCH ii regions (red star-like markers) and 10 in- 
frared protostars (red filled circles) have high luminosity. A few 
infrared-quiet cores (blue open circles) are weak MSX sources 
(like DR21(OH)) but most of them remain undetected (those 
more massive than 40 M© are plotted here). 



The Cygnus X molecular complex has already formed gen- 
erations of high-mass stars, since it contains several OB as- 
sociations ( Uyaniker et al. 2001)) and numerous H ii regions 
dWendker et al. 1991| l. We here discuss its ability to form high- 
mass stars in the near future by making a census of embed- 
ded H II regions (Sect. 15.2.11 ). high-luminosity protostellar cores 
(Sect. |5.2.2t . and massive infrared-quiet cores which are bright 
in SiO (Sect. |5.2.3t . We summarize in Table |5] the number of 
objects found at various stages of the high-mass star formation 
process. 



5.2.1. Embedded H ii regions 

Within our sample of high-density cores, we identify 15 
embedded H n regions from their bright infrared and cen- 
timeter continuum emissions (see Cols. 8-9 of Table [TJ. 
Most of them are already recognized as compact or 
ultra-compact H ii regions ( IDownes & Rinehart 19661 
IHaschick et al. 19811 IKurtz et al. 19941 IMolinari et al. 1998bl 
ITrinidad et al. 20031 ). Five ofliers (CygX-N47/DR21-D, 
CygX-N58/IRAS 20375-1-4109, CygX-S18/S106-IR, S20, 
CygX-S27/MSX 77.9550-h0.0058) coincide with radio cen- 
timeter sources from the 1.4 GHz Galactic Plane Survey 
dZoo nematkermani et al. 1990] ) or old centimeter maps 
dHarris 1973, Piphere tal. 1976[ ). As expected, observa- 
tions of high-density tracers such as CS, NH3 line and 350 yum 
continuum made toward these 1.2 mm sources suggest some 
free-free contamination (e.g. Motte et al. in prep.). 

Despite harboring an H 11 region, most of these massive 
young stellar objects are embedded and some of them (e.g. 
CygX-N46/DR21, CygX-N30AV75(N), CygX-S/AFGL 2591) 
even show high levels of ongoing star formation activity, such 
as maser emission and/or powerful outflow shocks (see Tables [T| 
and[3]and Figs.|5]and[8]). In agreement with these H 11 regions be- 
ing compact, ultra-compact, or even hyper-compact, our 1.2 mm 
study is sensitive to young and thus "embedded H 11 regions". 



5.2.2. "High-luminosity infrared sources" 

Dense cores which are luminous at infrared wavelengths are 
usually considered to be the best candidates to be (or to 
host) high-mass protostars (e.g. Wood & Churchwell 1989). 
Following a similar philosophy, we qualify as "high-luminosity 
infrared sources" those Cygnus X dense cores with bolomet- 
ric luminosity larger than 10^ Lq, which corresponds to that 
of a B3 star on the main sequence. This luminosity converts 
into a MSX 21 /vm flux of ~ 10 Jy, assuming that the lu- 
minosity of high-mass protostars is dominated by their mid- 
to far-infrared luminosity and their average colors are de- 
fined by [Wood & ChurchweU (1989)1 (see their Table 1). In 
practice, we estimate the total flux in the IRAS bands us- 
ing the simplified formulaeQ of jCasoli et al. (1986) and the 
|Wood & Churchwell (19 89) colors to express this total flux as 
a function of the single IRAS 25 fim flux. We finally convert the 
25 jjm flux into an expected MSX flux at 21 jjm by an inter- 
polation between 25 and 12 fim using the Wood & Churchwell 
color For d - 1.7 kpc, we finally estimate a typical MSX 
flux at 21 fim of ~ 10 Jy x (L/aAs/lO"' Lq) for a B3 protostar. 
Figure |7] displays the 21 yum flux detected toward Cygnus X 
dense cores as a function of their mass. With a 10 Jy limit in 
Fig. I2I we detect all the embedded H n regions plus 10 sources 
which are "high-luminosity infrared protostellar cores". Three 
of the latter are associated with OH, H2O, and/or CH3OH maser 
emission: CygX-N43AV75S-FIRl, CygX-N44/DR21(OH), and 
CygX-NW5A.kA225S. 

The high-luminosity infrared sources of the present 
Cygnus X survey are mostly IRAS sources fulfilling the 
Woo d & Churchwell (1989)1 criteria for UCH 11 regions or ob- 
jects in an earlier stage of high-mass star formation. Many of 
them were already identified as high-mass protostellar sources: 
e.g. W75S-FIR1 and -FIR2 (see [Chandler et al. 1993a| ) and 
IRAS 20343H-4129 and IRAS 20216-H4107 (Beuther et al. 2002). 



FjRAS - 1.75xlO-'3Wm-2x(5 12/0.79+525/2+560/3.9+5 100/9.9) 
with 5 12, 525, 560, 5 100 expressed in Jy. 
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Fig. 8. Integrated intensity of SiO(2-l) detected toward the 
Cygnus X dense cores more massive than 40 Mq, as a function of 
their mass. Infrared-quiet cores (blue open circles) are definitely 
stronger in SiO than the high-luminosity protostellar cores (red 
filled circles) and UCH n (red star-like markers). 



The profile of the SiO lines provides us with useful addi- 
tional constrains. Half of the infrared-quiet protostellar cores 
in Table [3] (CygX-N53, N63, N68, S30, S32 and possibly 
N3, N65, S7) show clear and smooth outflow wings in their 
SiO line profiles. These profiles are so clean that it can be 
taken as a good indication that the SiO emission is domi- 
nated by a single SiO outflow which, in turn, has to be pow- 
ered by a high-mass protostar. Four infrared-quiet protostellar 
cores (CygX-N44/DR21(OH), CygX-N48/DR21(OH)-S, CygX- 
N38/DR21(OH)-W, and CygX-N12) have a profile which is 
more indicative of two emission components: a strong narrow 
line plus a line with large wings (see Sect. 14.2b . This can arise 
from two different SiO outflows which would both still be pow- 
erful, or from a single outflow composed of two main veloc- 
ity patterns in the SiO shocks. An alternative interpretation is 
that, besides the outflow, some shocks associated with a hot core 
would be detected as a narrow SiO line component. 

Therefore, the > 40 Mq dense cores of Cygnus X all have a 
high probability of containing (at least) one high-mass protostar, 
whether or not they are bright infrared sources. Our selection 
of extreme-density cores has thus proven to be very efficient in 
locating sites of future formation of high-mass stars. 



5.2.3. Strong SiO emission as a probe for high-mass 
protostars 

In contrast to high-luminosity infrared sources, massive 
infrared-quiet cores are not all expected to host high-mass proto- 
stars. It is therefore interesting and somewhat surprising to find 
that all of them are bright SiO emitters (cf. Tableland Fig. |5]l. 
As shown below, their SiO(2-l) lines are well interpreted as due 
to a powerful outflow that is probably driven by one high-mass 
protostar. 

As shown in Fig. |8] which presents the SiO intensity of 
the Cygnus X dense cores versus their mass, the SiO lines de- 
tected for infrared-quiet cores are several times brighter than 
those observed for the high-luminosity sources. As an exam- 
ple, the well-recognized UCH n region CygX-S26/AFGL 2591 
is weaker (in SiO) than a completely unknown source such as 
CygX-N53. Since (SiO) outflows are tracing the accretion pro- 



cess (see Bontemps et al. 1996 1, this result suggests that massive 



infrared-quiet cores are in a more active accretion phase than 
high-luminosity infrared sources. Furthermore, the SiO emis- 
sion of massive infrared-quiet cores is several times stronger 
than the most extreme SiO shocks associated with nearby low- 
mass class protostars (see e.g. Table |3]l. This suggests that the 
infrared-quiet cores host either a clear single SiO outflow driven 
by a high-mass protostar, or a group of SiO outflows as extreme 
as those of L1448-MM or L1157-MM. Quantitatively the mean 
SiO brightness of the massive infrared-quiet cores in Cygnus X 
would require the summation of two low-mass class Os driving 
extreme SiO outflows (and up to 5 in the case of CygX-N53). 
This is rather improbable, since we evaluate the occurence of 
such SiO outflows to be of the order of 6% (2 out of 35 class Os 
known at d < 500 pc; e.g. lAndre et al. 2000] l. Statistically, the 
summation of several tens of low-mass class Os would be re- 
quired to explain the bright SiO emissions observed here. The 
SiO outflows of these infrared-quiet cores better resemble that of 



the high-mass protostar Orion-IRc2 (see Ziurys & Friberg 1987 
and Table O. These intensity comparisons suggest that the 
infrared-quiet protostellar cores of Cygnus X are much more 
active (in star formation) than nearby low-mass protostars, and 
that a large fraction of them should host high-mass analogs of 
class Os. 



5.3. "Massive infrared-quiet protostars": definition and 
lifetime 

5.3.1. Definition 

We define as "massive infrared-quiet protostellar cores" the 
Cygnus X dense cores that are more massive than 40 Mq, weak 
at infrared wavelengths (i.e. < 10 Jy at 21 fj.m), but present 
some signposts of stellar activity. Table [T] provides a list of 
33 dense cores more massive than 40 Mq (see also Table. |4| 
and Fig. [T] allows the identification of "high-luminosity in- 
frared sources" (see also Sect. I5.2.2] |. The 17 remaining dense 
cores should either be infrared-quiet protostellar cores or star- 
less cores. Among them, four are already identified as proto- 
stellar cores by their weak MSX detection and/or their associ- 
ation with masers: CygX-N44/DR21(OH), CygX-N65, CygX- 
S37/IRAS 20305-H4010, and CygX-N53. The majority of the 
17 massive infrared-quiet cores have been surveyed in SiO(2- 
1); 14 positions are listed in Table[3]and CygX-N32 is included 
in the beam targeting CygX-N30. The brightness and line shape 
of their SiO emission prove that the infrared-quiet cores are all 
protostellar in nature and that they probably host a high-mass 
protostar (see Sect. 15 .2. St . 



5.3.2. Lifetime of high-mass protostars 

Our study supports the existence of massive infrared-quiet pro- 
tostars, overlooked by infrared-based surveys (see Figs.|7}{8]l. It 
provides a sample of (mostly) new massive infrared-quiet pro- 
tostellar cores which contain high-mass protostars. We still lack 
the necessary angular resolution to identify individual high-mass 
protostars, but the characteristics of each Cygnus X dense core 
more massive than 40 Mq are believed to be shaped by those of 
a single embedded high-mass protostar (cf. Sects. l5n and l5.2.3] l. 
With this assumption, we give the first lifetime estimates for this 
new population of high-mass protostars, assuming a constant star 
formation rate in Cygnus X over the past 1 -2 Myr. 

Since our sample of Cygnus X dense cores is complete 
for embedded, massive, young stellar objects, we can statisti- 
cally estimate the relative lifetimes of massive infrared-quiet 
protostars and high-luminosity infrared sources (including in- 
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Table 5. Massive young stellar objects found in Cygnus X at various stages of the high-mass star formation process 





Starless Pre-stellar 
clumps cores 


Massive IR-quiet High-luminosity 
protostars IR protostars 


H II OB stars 
regions 


Number in Cygnus X 
References 
Statistical lifetime" 

< «H2 > 

Free-fall time* 


10 < 1 
this paper this paper 
~ 7 X 10' yr < 8 x 10^ yr 
6.7 X 10^^ cm-^^ ~ 1.5 X 10' cm"^^ 
4 X 10' yr ~ 9 x lO"* yr 


17 25 
this paper this paper 
-LSxlCyr -l^xW^yr 
1.7 X 10' cm-^^ 
8 X lO"* yr 


~ 800 2600 + 1000 
IWendker et al. 19911 IKnodlseder 20001 

~ 6 X 10' yr (2 ± 1) x 10*^ yr 


Number in Orion 


? ? 


^1 ^1 


~ 6 -81 



(a) Lifetime of massive starless clumps, pre-stellar cores, infrared-quiet protostars, high-luminosity infrared sources, and H ii regions measured 
relative to the known age of OB stars (cf. lHanson 2003b using the census given in Line 3. 

(b) Free-fall time measured from the mean values of the volume-averaged density given in Line 6: 'fj-ge-fall ~ -J 32 g'<p> ' 



frared protostars and embedded H ii regions). The separation 
of infrared sources between embedded H n regions and high- 
luminosity infrared cores is not final since centimeter continuum 
follow-up observations are necessary, and a low-level of cen- 
timeter free-free detection might be indicative of outflow shocks 
rather than H ii regions. In the following discussion, we therefore 
treat cores hosting high-luminosity infrared sources and possibly 
protostars as a single population. According to the numbers of 
high-luminosity infrared sources and massive infrared-quiet pro- 
tostars given in Table |5] (and taken from Sects. 15.2. H4sX2l and 
Sect. 15.3.1b . the infrared-quiet protostellar phase of high-mass 
stars might last as long as their better-known infrared-bright pro- 
tostellar (or embedded H n) phase. This result is in disagreement 
with the idea that high-mass protostars are very luminous from 
the very first stage of their evolution and thus that their infrared- 
quiet phase is very short-lived. 

Table |5] gives our current estimates of the number of mas- 
sive young stellar objects forming and stars already formed in 
Cygnus X. The number and age of OB stars in Cygnus X are 
taken to be those of the Cyg OB2 association (IKnodlseder 2()()()i 
IHanson 20031 1. increasing uncertainties to take into account 
the other (smaller) OB associations which have different ages 



(Uyaniker et al. 2001 1. The number of H ii regions is crudely 



given by the numbers of centimeter free-free sources observed 
in Cygnus X (Wendke r et al. 199l1 l. The statistical lifetimes of 
high-mass protostars derived relative to the number of OB stars 
give a time range of (0.5 - 3.1) x 10"^ yr for the infrared-quiet 
protostellar stage and (1.2 - 7.9) x 10^ yr for the complete pro- 
tostellar stage (see also Table|5]l. For the sake of comparison, the 
statistical lifetime measured for the high-mass protostellar phase 
in Cygnus X (~ 3.2 x 10"* yr) is one order of magnitude smaller 
than the typical age of nearby, low-mass class I protostars (e.g. 



Kenyon & Hartmann 1995 1 and more closely resembles that of 
low-mass class Os (cf. I Andre et al. 2000] l. 

Figure |9] displays the number of high-mass protostars, H ii 
regions and OB stars (i.e. objects in the three main evolutionary 
phases) as a function of the time elapsed since the beginning of 
the protostellar collapse (i.e. since the formation of a hydrostatic 
stellar embryo). Absolute lifetimes of the high-mass protostel- 
lar phase (encompassing both the infrared-quiet and the high- 
luminosity infrared protostellar phases) as well as those of the 
high-mass starless phases can be roughly estimated from their 
free-fall dynamical timescales. Table |5] gives the estimates we 
made using the mean values of the volume-averaged density of 
massive dense cores observed in Cygnus X. The most massive 
dense cores of Cygnus X should thus collapse in a free-fall time 
of ~ 8 X lO'* yr, which is taken to be the end of the high-mass 
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Fig. 9. Number of massive young stellar objects as a function 
of the time elapsed since the beginning of the protostellar col- 
lapse in Cygnus X and Orion. The age estimates (also given in 
Sect. 15.3.2b are not statistical, in order to have the X-axis and 
Y-axis independent. 



protostellar phase in Fig.|9] The H n regions phase should last for 
several times the typical lifetime of UCH ii regions (~ 10^ yr, cf. 
Wood & Churchwell 1989), starting in Fig. ||]at the end of the 
protostellar phase. The free-fall time of high-mass protostars in 
Cygnus X agrees, within a factor of 3, with the statistical life- 
times of the complete population of high-mass protostars (in- 
cluding both infrared-quiet and high-luminosity infrared proto- 
stars) derived relative to the number of OB stars (see Table|5]l. 

To compare Cygnus X to Orion, Table |5] makes the cen- 
sus of all massive young stellar objects known in Orion, and 
Fig. |9] shows their distribution as a function of time. The num- 
ber of high-mass stars (down to B3 spectral type) within all 
the Orion OB associations is measured by |Brown et al. (1994)] 
while the H ii regions are associated with the well-known reflec- 
tion nebulae M42, NGC 2024, NGC 1977, M43, NGC 2068, and 
NGC 2071. One source could qualify as a high-luminosity in- 
frared protostar (Orion-IRC2) and another as a massive infrared- 
quiet protostar (OMCl-S). Interestingly, the massive young stel- 
lar objects forming, and stars formed, in Cygnus X are 20 times 
more numerous than in Orion (see Table |5]), in good agreement 
with the mass ratio (13) of the Cygnus X and Orion cloud com- 
plexes. 
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5.4. Where is the pre-stellar phase of massive star 
formation? 

Among the 33 most massive (> 40 Mq) dense cores of 
Cygnus X, no starless cores are found, with the possible ex- 
ception of CygX-N69 which was not selected for SiO observa- 
tions because of its relatively large size and corresponding low- 
density (~ 10"* cm~^). This non-detection of starless cores means 
that high-mass pre-stellar cores, if they exist, are short-lived 
(< 10^ yr, see Table|5]l. This is particularly surprising since the 
lifetime of pre-stellar cores and condensations in nearby, low- 
mass star-forming regions is more than two orders of magnitude 
longer (1 - 4 x 10 ^ yr, e.g. lOnishi et al. 2002t iKirk et al. 20051 
I Andre et al. 2007l l. Below, we consider two interpretations: ei- 
ther the star formation activity of Cygnus X is not continuous 
(see Sect. 15.4. Il l or massive pre-stellar clumps are dynamically 
evolving into protostars (see Sect. 15.4.21 ). 



5.4.1 . Star formation in Cygnus X cannot be caused by a 
burst 

If the Cygnus X molecular complex were experiencing a sudden 
burst of star formation, it might be due to the interaction with 
the H II regions and/or OB associations, including the Cyg OB2 
cluster In this case, we expect differences in the concentration of 
gas and in the population of young stellar objects in the CygX- 
North and CygX-South regions. 

Evidence shows that the concentration of gas is differ- 
ent in CygX-North and CygX-South. The CO survey of 
[Schneider et al. (2006)1 the present MAMBO-2 maps, and 
our small-scale analysis trace various cloud structures (see 
Sects. [m43. 21 ). We can thus measure the gradual concentration 
of molecular clouds (0.5-50 pc scales), first into 0. 1-5 pc scales, 
and then into 0. 1 pc dense cores. In CygX-North (280 x 10"* Mq, 
cf. Schneider et al. 2006), the concentration factor is 8 for each 
step while in CygX-South (480 x 10^ Mq) it is 18 for each step. 
Therefore, when considering the mass involved in dense cores, 
the CygX-North region is 5 times more concentrated (CO clouds 
concentrate into 0.1 pc cores with an efficiency of ~ 1.5%) than 
CygX-South (concentration efficiency of ~ 0.3%). This suggests 
that CygX-North could be 5 times more efficient in forming 
dense cores, and possibly forming stars, than CygX-South. 

In contrast, the relative numbers of "high-luminosity infrared 
sources" versus "massive infrared-quiet protostellar cores" are 
identical (within our large error bars) in both regions. Moreover, 
other generations of (massive) young stars have already formed 
nearby (see e.g. the MSX sources at ~ 0.2 pc from the dense 
cores). The agreement between the free-fall lifetime and the sta- 
tistical lifetime of high-mass protostars relative to OB stars also 
suggests that star formation has been continuous, at least over 
the past ~ 10^ yr (see also Fig.|9]). Therefore, it is highly im- 
probable that Cygnus X is currently experiencing a burst of star 
formation. 



5.4.2. Is tine star formation in Cygnus X dynamic? 

The present study fails to find a single good candidate for be- 
ing a pre-stellar dense core, i.e. a starless cloud structure with 
a ~ 0.1 pc size and a ~ 10^ cm"^ volume-averaged density, 
which is gravitationally bound. This does not preclude, how- 
ever, the existence of smaller and denser pre-stellar condensa- 
tions (~ 0.01 pc size and ~ 10^ cm"-' density according to Motte 
et al. 1998, see Tabled within the Cygnus X protostellar cores. 
On the other hand, we have identified 10 clumps which may be 



starless (~ 0.8 pc size and 7 x 10^ cm"^ density, see Tables |2] 
and nil. The Virial mass of these clumps (Myir - 3Rcr^/G, 
Bertoldi & McKee 1992 ) are estimated from the NaH^ linewidth 
measured at their peak location (Schneider et al. in prep.) and 
the FWHM size given in Table |2l they suggest that most of these 
clumps are close to being gravitationally bound. 

From the numbers of starless clumps and pre-stellar cores 
detected in Cygnus X and compared to the numbers of OB 
stars, we statistically estimate their lifetimes to be ~ 10"^ yr 
and < 10^ yr, respectively (see Table |5]l. These lifetime values 
are considerably smaller (by factors of 60-110) than the free- 
fall times we estimated. They are also shorter than the lifetimes 
observed in nearby regions that form mostly low-mass stars, 
since starless cloud structures with volume-averaged densities 
of ~ lO"' cm"^ and ~ 10'* cm""* have pre-stellar lifetimes of 
~ 10^ yr and ~ 10^ yr (see Fig. 11 of IKirk et al. 2005] and ref- 
erences therein). Therefore, massive pre-stellar cores seem to be 
short-hved or even transient cloud structures compared to nearby 
low-mass pre-stellar cores, to high-mass protostars, and to mas- 
sive starless clumps. 

A dynamical process should thus be acting to create pre- 
stellar condensations from starless clumps. It should be super- 
sonic to concentrate mass from densities of 7 x lO-' cm"-' to 
2 x 10^ cm"^ in a time as short as, or even shorter than, the 
free-fall timescale. Similarly short lifetimes are theoretically ex- 
pected in molecular clouds where highly turbulent processes 
dominate and pre-stellar cores are magnetically supercritical 
(e.g. |Vazquez-Semadeni et al. 2005 1. Such dynamical processes 
are also necessary for the protostellar lifetime to last for only one 
free-fall time (see Table |5]l and, consequently, for the accretion 
to be strong enough (~ 10"^ Mq yr"') to overcome the radia- 
tion pressure and form a high-mass star. In fact, the DR21 fila- 
ment is observed to undergo a global supersonic coUapse, possi- 
bly due to a strong external compression (IChandler et al. 1993bt 
Mot te et al. 20051 ). This is encouraging but it remains to be veri- 
fied whether such a dynamical process acts throughout the com- 
plex. 

Taken together, the lifetimes measured here suggest a faster 
evolution (by a factor 10 at least) of the pre-stellar and proto- 
stellar phases of high-mass stars compared to those typically 
found for nearby, low-mass stars. Molecular clouds harboring 
high-mass star formations should thus be quantitatively different 
in terms of e.g. density and kinematics than those of our close 
neighborhood. 



6. Conclusions 

In order to improve our knowledge of the high-mass star for- 
mation process, we started an unbiased search for its earliest 
phases, i.e. the high-mass analogs of low-mass pre-stellar cores 
and class protostars. One of the best targets for such a study is 
Cygnus X, which we recently recognized as the richest molec- 
ular and H ii complex at less than 3 kpc (it is located at 1 .7 kpc 
from the Sun). We have therefore made a complete 1.2 mm con- 
tinnum survey of the Cygnus X cloud complex and have per- 
formed SiO(2-l) follow-up observations of the best candidate 
progenitors of high-mass stars. Our main findings can be sum- 
marized as follows: 

1. Our MAMBO-2 (1.2 mm, 11" resolution) imaging of the 
Cygnus X molecular cloud complex gives a complete view 
of the cloud structures ranging from 0.03 to 5 pc, i.e. from 
dense cores to clumps (cf. the terminology used for nearby 
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molecular clouds). We performed a multi-resolution analy- 
sis to identify a few tens of massive large-scale (~ 0.7 pc) 
clumps and extract 129 compact (~ 0.1 pc) dense cores. 

2. The Cygnus X dense cores are the high-density parts of mas- 
sive clumps which resemble HMPOs or IRDCs. The dense 
cores have similar sizes (~ 0.1 pc) but mass and volume- 
averaged densities (4 - 950 Mq, ~ 1.1 x 10^ cm"-*) larger 
than those of nearby, pre-stellar dense cores. The most mas- 
sive (> 40 Mq) and thus "extreme-density" cores of this large 
and unbiased sample are probable precursors of high-mass 
stars spanning, for the first time, the whole evolution from 
pre-stellar cores to "embedded H ii regions". 

3. We have used the MSX 21 jjm flux arising from our 
MAMBO-2 dense cores to identify the high-luminosity (> 
lO-' Lo) massive young stellar objects of Cygnus X: 15 
UCH II regions and 10 "high-luminosity infrared protostars". 
Strikingly, half of the dense cores considered as the best can- 
didate precursors of high-mass stars are found to be infrared- 
quiet (i.e. 17 Cygnus X dense cores more massive than 
40 Mo have f 2i^m < 10 Jy). 

4. We have surveyed the most massive (> 40 Mq) dense cores 
of Cygnus X in SiO(2-l) to search for shocked gas in molec- 
ular outflows and/or hot cores. The association of high- 
velocity SiO emission with all massive infrared-quiet cores 
provides persuasive evidence that stars are already forming 
in these cores as outflows are tracing the accretion process. 
The brightness and shape of the SiO line compared to that 
of nearby, low-mass protostars suggest that these extreme- 
density cores may host "massive infrared-quiet protostars". 

5. Our unbiased survey of the massive young stellar objects in 
Cygnus X shows that massive infrared-quiet protostars do 
exist, and that their lifetime should be comparable to that of 
more evolved high-luminosity infrared protostars. By com- 
paring the number of high-mass protostars and OB stars in 
the whole of Cygnus X, we estimate a statistical lifetime of 
3 X 10"^ yr for high-mass protostars. One order of magni- 
tude smaller than the lifetime of nearby, low-mass protostars, 
such a value agrees, within a factor of 3, with the free-fall 
time of Cygnus X dense cores. 

6. Our complete census of massive young stellar objects in 
Cygnus X fails to discover the high-mass analogs of pre- 
stellar dense cores (~ 0.1 pc, > lO'* cm"^). Their correspond- 
ing lifetime (< 10^ yr) is smaller than one free-fall time, 
in marked contrast with nearby, pre-stellar cores that typi- 
cally live for 2-5 free-fall times. We propose that the star- 
less, massive but lower-density clumps (800 Mq, ~ 0.8 pc, 
~ 7 X 10^ cm""*) that we observe rapidly concentrate and 
collapse to form high-mass protostars. 

7. Lifetime measurements of the pre-stellar and protostellar 
phases of high-mass stars in Cygnus X suggest a dynam- 
ical contraction for pre-stellar cores and a supersonic evo- 
lution for protostars. Highly turbulent processes throughout 
the Cygnus X molecular cloud complex would be necessary 
in such a dynamical picture of the high-mass star formation 
process. 

8. Our Cygnus X study suggests that far-infrared to sub- 
millimeter continuum imaging of entire star-forming com- 
plexes (such as those proposed with Herschel by Motte, 
Zavagno, Bontemps et al.: the HOBY^ Key Programme) 



' "HOBYS: the Herschel imaging survey of OB Young Stellar ob- 
jects" is a Guaranteed Time Key Programme jointly proposed by the 
SPIRE and PACS consortia, and the Herschel Science Centre. 



will definitely contribute to a better knowledge of the high- 
mass star formation process during its earliest phases. 
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Appendix A: A skynoise reduction technique dedi- 
cated to extended sources 

Ground-based (sub)millimeter continuum observations of astro- 
nomical sources are hampered by the emission of the atmosphere 
and its fluctuations, called "skynoise". The bulk of the atmo- 
spheric emission and its slow variations in time can be removed 
by the dual-beam scanning mode generally used for bolome- 
ter observations (the "EKH method", lEmerson et al. 1979] !. 
However, the residual atmospheric fluctuations generate an ex- 
cess noise in the restored maps that can dominate the Gaussian 
noise (see e.g. Fig.fTOk). In the present study, more than half of 
the on-the-fly maps were observed with a skynoise level quali- 
fied as medium or high. In our maps, medium skynoise implies 
cr - 30-60 mJy beam"' and high skynoise cr > 60mJybeam"', 
while the nominal rms is cr ~ 15 mJybeam"'. This excess 
noise can theoretically be suppressed by using the fact that it 
is well correlated across the whole array. Several skynoise re- 
duction algorithms have been developed assuming that either 
the atmospheric emission is equal in all receivers at all times 
("mean" or "neighbor" approximation), or that it is correlated 
on a few seconds of time with correlation coefficients decreas- 
ing when the angular separation of receivers increases (see e.g. 
NIC and MOPSIC IRAM softwares for bolometer-array data, 
Broguiere et al. 1995| and |Zylka 1998| . While efiicient for the 
detection of pointlike sources, such assumptions can have dra- 
matic consequences on the representation of the extended emis- 
sion. Indeed, the large-scale structure of a source leads to a corre- 
lation of bolometer signals which is non-zero, and caution must 
be taken not to filter out most of the extended emission. 

We have developed a program which uses a source model 
and the mean approximation to remove the skynoise while pre- 
serving, as much as possible, the structure of sources like molec- 
ular clouds (compare Figs. [TOb-b). We describe below its four 
steps: 

(1) The first step consists of creating a good source model. A 
first-order data processing is done by reducing all individual 
maps and combining them with weights consistent with their 
intrinsic quality, i.e. accounting for both Gaussian noise and 
skynoise levels (see e.g. Fig. [TOb). This first-order mosaic is 
used as a model of the source after cutting its noisier parts 
(mostly map edges) and removing artefacts such as negative lev- 
els stronger than the rms. 

(2) This source model is then subtracted from all the original 
scans. To do so, we first simulate the on-the-fly dual-beam sig- 
nals of each bolometer when the MAMBO-2 camera is mapping 
the source model with the observational geometry of the real 
scans (see IMotte & Andre 20011 and the routine "simulate" in 
NIC). These simulated source signals are then removed from the 
original on-the-fly dual-beam signals recorded by each bolome- 
ter. At this point, and for each scan, the bolometer data streams 
should only contain Gaussian noise and skynoise. 

(3) As a third step, the skynoise emission contained in each 
processed scan is estimated by using the averaging method. 
While cruder than any "correlation method" (like that proposed 
in NIC), it allows the treatment of each observing point (taken 
every half a second) independently. We often find the averaging 
method to be more eflicient, suggesting that the skynoise power 
is strong at high frequency. 

(4) Afterwards, the original scans are reduced for skynoise by 
subtracting, for each scan, the skynoise data stream estimated 
above from the on-the-fly dual-beam signals of each bolometer. 
Finally, the individual scans, now treated for skynoise, are re- 




s;^ Before skynoise reduction 




Ml 

7^ After skynoise reduction 



a (2000) 

Fig. 10. 1.2 mm map of the DR13 filament in the Cygnus X 
molecular cloud complex obtained with MAMBO-2 at the 
IRAM 30 m telescope. Contour levels are -45 mJy beam"' and 
then 45 to 315 mJybeam"' with steps of 45 mJybeam"'. 
(a) Mosaic built from individual maps taken with low to 
medium skynoise. As a result, the rms noise level ranges from 
15 mJybeam"' to 60 mJybeam"'. (b) Mosaic built from in- 
dividual maps processed by our skynoise reduction program. 
The resulting noise level is relatively uniform with cr - (14 + 
2) mJy beam"'. 



duced with the default procedure and combined to create a better 
(second-order processed) mosaic (see e.g. Fig.fTOb). 

This process can evidently be iterated to improve the map 
sensitivity but a single iteration already reduces most of the 
skynoise and probably gives more secure results. The final 
individual maps have a noise rms of 20 + 5 mJybeam"', 
15 mJybeam"' being the nominal rms in the absence of 
skynoise, and with a zenith atmospheric opacity of ~ 0.2. Our 
skynoise reduction technique has thus improved the map sensi- 
tivity by a factor of ~ 2 on average, and up to a factor of 6 for 
the poorer quality original maps. 

Our method relies on the quaUty of the source model and 
thus mainly on the redundancy of observations taken in every 
single pixel. The redundancy is currently ~ 20 for one single 
map and ~ 40 on average, since most parts were imaged sev- 
eral times. The mosaics shown in Figs. |2b-c recover most of the 
extended emission up to ~ 10', which is the azimuthal size of 
individual maps. In agreement, the Cygnus X main region has 
recently been imaged with the SHARC-II camera installed at the 
CSO telescope (Motte et al. in prep.). Those data generally con- 
firm the cloud structures detected in our MAMBO-2 maps. We 
acknowledge, however, that the extended emission is attenuated 
in a couple of sites where data are too noisy and not redundant 
enough to build a good source model (e.g. south of DR23). 
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Appendix B: An extraction technique dedicated to 
compact sources 

Since stars are generally forming in the densest parts of molec- 
ular clouds, we have developed a source extraction technique 
aiming at identifying and extracting dense compact fragmeri 
that should be the best potential sites of star formation. 

In spirit, this method is equivalent to an eye search of 1 
cal density peaks followed by flux measurements with an ape 
ture optimized for each source. We outline below the procedu 
for this automated method, which was originally developed 1 
Motte et al. (2003): 

( 1 ) As a first step, one has to determine the cutoff" lengthsca 
at which the sources of interest (here dense cores, cf. Fig. [TT| 
separate from their surrounding (here clumps and clouds, c 
Fig. [TTb). This scale corresponds to twice the mean outer r 
dius of the most compact objects one can distinguish on tl 
MAMBO-2 map or the mean projected distance between thei 
This scale is typically about 5 to 10 times the HPBW. 

(2) For further analysis, we only keep spatial scales smaller th; 
the cutoff" lengthscale determined at step 1 (here 1 pc). For this 
purpose, we use a multiresolution program based on wavelet 
transforms (MRE, cf. Starck & Murtagh 2006). This makes a se- 
ries of filtering operations and provides "views" of the image 
at different spatial scales. We therefore surr0 all wavelet planes 
of the original image up to that scale to create images like that 
shown in Fig. [TTb . 

(3) As a third step, we use a version of the Gaussclumps pro- 
gram dStutzki & Glisten 19901 IKrameret al. 19981 which has 
been customized for continuum images (see I Motte et al. 20031 
and Mookerje a et al. 2004[ ) and optimized for compact source^. 
Gaussclumps is applied to the filtered image created at step 2. 
It identifies compact fragments above a 5 cr level (here 
75 mJybeam"') and gives first-guess parameters of a 2D- 
Gaussian fit. If the Gaussian shape is not ideally fitting com- 
pact fragments, it gives appropriate estimates of their size and 
mass even when their environment is not perfectly filtered 
out. Gaussclumps is known to be efficient at deblending com- 
pact sources. Noisy edges that usually pose problems (see e.g. 
IRrameretal. 1998 1) have been removed. 

(4) Finally, the most critical and time consuming step is to im- 
prove the extraction of individual fragments from their surround- 
ing background. The main goal is to avoid Gaussian fits being 
biased to a single lengthscale, which would be that determined 
at step 1 . We therefore repeat steps 2 and 3 for the sources iden- 
tified at step 3 and use a cutoff scale optimized for each frag- 
ment. We estimate the scale up to which one fragment is resolved 
from its surrounding to be twice the first-guess measure of its 
Gaussian FWHM (see step 3). Therefore, for each of the (here 
129) sources, we create a final-guess filtered image and adjust 
its Gaussian fit parameters. 




Fig. 11. Separation of cloud structures from the MAMBO-2 map 
of S106 shown in (a) into (b) dense cores (i.e. compact struc- 
tures with size scales smaller than 128") and (c) their environ- 
ment (i.e. cloud structure with a typical scale > 128"). Dense 
cores extracted by the technique described in Appendix B are 
indicated by crosses. 



The original image can be expressed as the sum of all wavelet 
views (also called planes) plus the smoothed image (last plane contain- 
ing all remaining scales). 

' Gaussclumps parameters for initial guesses are: FWHM= 1.1 x 
HPBW, aperture = 1.5 x HPBW 
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Table 1. Properties of dense cores (compact cloud fragments) detected in the Cygnus X complex 





Fragment 


^ peak a 
1.2 mm 




^ int a 
1.2 mm 




< >■= 


MSX'' 


Name, cormnents, and coincidence 




name 


(mJy/ 
beam) 


(pc X pc) 


(mJy) 


(Mo) 


(cm-3) 


21 fim 

(Jy) 


with stellar activity signatures 


CygX-North region 


Nl 


J203532 7+421953 


120 


15 X 09 


320 


17 


4.7 X 10* 






N2 


J203534.0+422023 


110 


0.22 X 0.07 


380 


20 


4.0 X lO'* 


3.74 


MSX 81.3039+1.0520 


N3 


J203534.2+422005 


720 


1 1 X 09 


1580 


84 


3.6 X 10' 






N4 


J203558. 0+422431 


90 


0.06 X 0.06 


120 


< 6 


< 1.2 X 10' 




probably not an embedded YSO^ 


N5 


J203606 5+413958 


160 


0.12 X 0.07 


350 


19 


9.4x10^* 






N6 


1203608 1+413958 


260 


23 X 1 4 


1280 


44-68 


3 2-5 OxIO'* 


90.88 


IRAS 20343+4129 


N7 


J203615 7+425305 


80 


0.12 X 0.10 


190 


10 


3.2 X 10* 






N8 


J203638.4+422911 


170 


17 X 09 


500 


27 


6.4 X 10* 


Abs 




N9 


J203639.5+425113 


100 


0.23 X 0.14 


500 


26 


1.7 X 10* 


Abs 




NIO 


J203652.2+413623 


590 


0.14 X 0.09 


1530 


<52-81 


<1. 5-2.4x10' 


99.33 


IRAS 20350+4126, cm' free-free, 
maser^ (CH3 OH/OH) 


Nil 


J203656.8+421320 


180 


0.11 X 0.06 


340 


18 


1.4 X 10' 


Abs 




N12 


J203657.4+421 127 


740 


0.11 x0.09 


1630 


86 


3.8 X 10' 


Abs 




N13 


J203658.4+421118 


90 


0.04 X 0.03 


100 


5 


7.1 X 10' 


Abs 




N14 


J203700.9+413457 


500 


10 X 07 


940 


50 


3.5 X 10' 


19.71 


IRAS 20352+4124 


N15 


J203710 9+413339 


80 


09 X 03 


110 


6 


2.1 X 10' 






N16 


J203716. 9+421631 


210 


0.12 X 0.12 


550 


29 


7.5 X 10* 


3.87 


MSX 81.4452+0.7635 


N17 


J203726 5+413319 


100 


0.15 X 0.03 


200 


10 


1.6 X 10' 






NIB 


J203730.2+421359 


110 


0.17x0.13 


430 


23 


2.7 X 10* 


- 




N19 


J203745. 8+424343 


80 


0.09 X 0.03 


110 


6 


1.9 X 10' 






N20 


J203747 .4+423842 


100 


10 X 03 


150 


8 


2.1 X 10' 


2.06 


MSX 81 7947+0 9116 


N21 


170^801 7+4?^9^S 


100 




380 


20 


2.5 X 10* 






N22 




90 


27 X 1 1 


420 


22 


1.9 X 10* 






N23 


T203804 Q+4233 1 7 


80 


22 X 06 


250 


13 


4.0 X 10* 






N24 


1203805 8+423952 


140 


17x017 


650 


34 


3.0 X 10* 






N25 


1203810 4+423805 


170 


12x011 


450 


24 


6.5 X 10* 






N26 


J203821.2+421131 


100 


0.20 X 0.10 


340 


18 


2.6 X 10* 






N27 


J20383 1.6+423804 


90 


0.10 X 0.10 


200 


11 


4.1 X 10* 






N28 


J203833 l+M'i^\A 


140 


0.19 X 0.11 


510 


27 


3.9 X 10* 






N29 


J203834.6+420619 


110 


0.25 X 0.13 


540 


29 


2.0 X 10* 






N30 


J203836.6+423730 


3650 


0.15 X 0.10 


10640 


<361-563 


<0.8-1.2xl0* 


136 


W75N(B), cm free-free^, 
maser* (CH3OH/H2O/OH) 


N31 


J203836.9+423750 


120 


0.03 X 0.03 


130 


7 


1.2 X 10** 




close to H n region W75N(A) 


N32 


J203838.3+423733 


540 


0.16x0.10 


1630 


86 


1.7 X 10' 


- 




N33 


J203847.0+423800 


100 


0.08 X 0.07 


170 


9 


8.1 X 10* 






N34 


J203854.3+421917 


150 


0.19 X 0.07 


430 


23 


6.1 X 10* 






N35 


J203854.4+421855 


90 


0.13x0.03 


160 


8 


1.5 X 10' 


- 




N36 


J203858.5+422330 


90 


0.19 X 0.03 


210 


11 


1.2 X 10' 






N37 


J203858.6+422435 


80 


0.03 X 0.03 


80 


4 


8.7 X 10' 






N38 


J203859.4+422223 


830 


0.18 X 0.09 


2610 


138 


2.7 X 10' 




DR21(OH)-W 


N39 


J203859.9+422732 


70 


0.16 X 0.08 


200 


10 


3.0 X 10* 






N40 


J203859.8+422343 


400 


0.33 X 0.08 


2000 


106 


1.1 X 10' 




DR21(OH)-N2 


N41 


J203900.2+422325 


180 


0.03 X 0.03 


180 


10 


1.9 X 10*" 




DR21(OH)-Nl 


N42 


J203900.3+421907 


90 


0.05 X 0.03 


100 


5 


4.5 X 10' 






N43 


J203900.6+422435 


370 


0.19 X 0.10 


1250 


66 


1.1 X 10' 


15.9 


W75S-FIR1, maser''* (H2O/OH) 


N44 


J20390 1.0+422246 


2970 


0.14 X 0.11 


8430 


446 


1.0 X 10* 


6.78 


DR21(0H), 

maser^'* (CH3OH/H2O/OH) 


N45 


J203901.1+421814 


130 


0.13 X 0.08 


300 


16 


6.4 X 10* 


Abs 




N46 


J203901. 4+421934 


4200 


0.19x0.14 


17930 


<609-949 


<5.9-9.1xl0' 


272 


DR21, cm'' free-free, 
maser« (HjO) 


N47 


J20390 1.4+42 1954 


630 


0.08 X 0.03 


840 


<28-44 


<1.2-1.9xl0* 


272 


DR21-D, cm'' free-free 


N48 


J20390 1. 5+422203 


1160 


0.17x0.11 


3720 


197 


3.5 X 10' 




DR21(OH)-S 


N49 


J203902.0+422700 


80 


0.15 X 0.03 


160 


8 


1.2 X 10' 






N50 


J20390 1.9+42 1835 


80 


0.24 X 0.03 


240 


12 


1.0 X 10' 






N51 


J203902.4+422459 


680 


0.19x0.10 


2370 


125 


2.0 X 10' 


26.50 


W75S-FIR2 


N52 


J203903.0+422615 


180 


0.16x0.12 


630 


33 


4.9 X 10* 






N53 


J203903.2+422549 


590 


0.14x0.09 


1600 


85 


2.2 X 10' 




maser^ (CHjOH/HjO) 


N54 


J203903.6+422530 


150 


0.18x0.08 


450 


24 


6.0 X 10* 


10.85 


W75S-FIR3 


N55 


J203906.1+421819 


100 


0.15 X 0.05 


220 


11 


8.0 X 10* 


Abs 




N56 


J203916.9+421607 


270 


0.12x0.08 


600 


32 


1.4 X 10' 


3.68 


MSX 81.6632+0.4651 


N57 


J203919.3+421556 


130 


0.09 X 0.05 


210 


11 


1.5 X 10' 


Abs 


EROl 
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Table 1. Continued 





Fragment 


g peak a 
1.2 mm 


FWHM^ 


5/"t " 
1.2 mm 


^1.2 mm' 


< KHz >' 


MSX'^ 


Name, comments, and coincidence 




name 


(mJy/ 
beam) 


(pc X pc) 


(nUy) 




(cm ) 


21 /jm 


with stellar activity signatures 


JN JO 






U.ZJ X U. IV 


T7^A 
Z/JU 


^no 1/1/1 


<4.z-o.dX1U 


o34. 13 


IRAb zU375+41i)y, cm iree-rree 


iNjy 


JZUjy J i .Z+4iZUUj 


ZUU 


U.lz X U.Uj 


Q7A 


on 
zU 


z.U X lU 






JNOU 


JzUjyjo.z-l-41 Lyil 




U.lO X U.l J 


^AA 


'id 

JU 


4.Z X ilT 






MAI 


TOniOCT ^lA 1 <oi 1 
JZiJjyJ / .o+^LjyLL 


lUU 


A 1 A AO 

U.lz X u.uy 


01A 
ZjU 


iZ 


4. / X llT 


ll.OJ 


ivioA. oi.jioo+u.iyzo 


INOZ 


JZU4UUJ.O+41Z/JO 


14U 


A 1 Q v A 1 
U. lo X U.IZ 


<AA 
jUU 


LI 


J.O X iir 




n^nna *n A jcv 01 11 AO A 1 /ten 

Close to M5>x isi.iiuy-u.i4Dy 

nr^A AjTCV 01 1 1OC A 1 1A1 

anaMbA 5l.lz2D-u.l343 




J ZU4UU J .Z+4 ijZlJ 


/4U 


n HQ V A A^ 

U.Ub X U.Uj 


1 AOA 

luyu 


JO 


1.1 X lU 


ADS 






J ZU4UZ / .Zi-4 i jOJ i 


1 nn 

lUU 


n 1 o V A AO 
u. iz X u.uy 


zzu 


IZ 


4.0 X lU 


/\DS 




N65 


J204028.4+415711 


230 


0.16 X 0.12 


780 


41 


6.4 X lO'' 


2.0 




N66 


1204029.8+414933 


110 


0.15x0.06 


250 


13 


6.1 X 10^ 


Abs 




N67 


J20403 1.0+414957 


100 


0.08 X 0.05 


150 


8 


1.4 X 10' 


Abs 




N68 


J204033.5+415903 


390 


0.10x0.10 


840 


44 


2.0 X 10' 


Abs 




N69 


J204033.7+415059 


110 


0.47x0.18 


1280 


68 


1.2 X lO'' 






N70 


J204034.3+413845 


80 


0.12x0.07 


170 


9 


5.0 X 10^ 


Abs 




N71 


J204038.2+4 15259 


80 


0.21 x 0.06 


230 


12 


3.4 X lO'' 






N72 


J204045.7+415750 


100 


0.28x0.13 


540 


29 


1.7 X lO'' 






CygX-South region 


SI 


J201643.3+392313 


90 


0.08 X 0.05 


140 


8 


1.1 X 10" 


3.08 


IRAS 20149+3913 


S2 


J201648.5+392212 


100 


0.17x0.06 


260 


14 


5.4 X lO'' 


Abs 




S3 


J201658.6+392106 


210 


0.12x0.08 


440 


23 


1.2 X 10' 


5.27 


IRAS 20151+3911 


S4 


J201737.0+392139 


90 


0.14x0.10 


230 


12 


3.1 X lO'' 


Abs 




S5 


J201745.3+392041 


120 


0.16x0.09 


340 


18 


4.1 X 10^ 






S6 


J202036.5+393821 


120 


0.03 X 0.03 


120 


6 


1.2 X 10* 






S7 


J202038.0+393820 


280 


0.22 X 0.03 


760 


40 


3.5 X 10' 






S8 


J202038.8+393751 


900 


0.19x0.12 


3540 


<121-188 


<1.4-2.1xl0' 


145.48 


Moll 2 IS, cm'" free-free, 
maser" (H2O) 


S9 


J202039.5+393738 


120 


0.03 X 0.03 


130 


7 


1.1 X 10*^ 






SIO 


J202044.4+393520 


180 


0.19x0.14 


760 


40 


4.1 X lO'' 






Sll 


J202153.7+395933 


110 


0.15x0.11 


340 


18 


3.6 X lO'' 






S12 


J202220.1+395821 


130 


0.16x0.08 


350 


18 


5.5 X 10* 




close to IRAS 20205+3948 


S13 


J202632.4+395721 


160 


0.15x0.11 


500 


27 


5.0 X 10* 


4.84 


MSX 78.3762+1.0191 


S14 


J202634.8+395725 


200 


0.09 x 0.05 


330 


17 


2.4 X 10' 


Abs 




S15 


J202714.0+372258 


120 


0.29x0.12 


630 


34 


2.2 X lO"* 


Abs 




S16 


J202717.2+372218 


260 


0.11 X 0.05 


460 


24 


3.0 X 10' 






S17 


J2027 18.7+372227 


130 


0.07 X 0.06 


210 


11 


1.4 X 10' 


_ 




S18 


J202726.0+ 372246 


960 


0.21 X 0.13 


4270 


< 145-226 


<1.2-2.0xl0' 


413.63 


S106-IR, cm'2 free-free, 
maser* (H2O) 


S19 


J202726.9+ 372233 


310 


0.08 X 0.03 


420 


< 14-22 


<6.1-9.5xl0' 


413.63 


part of S106-IR, cm'^ free-free 


S20 


J202727.4+372256 


510 


0.11 xO.ll 


1300 


<44-69 


< 1.2-2.0x10' 


413.63 


part of S106-IR, cm' free-free 


S21 


J202729.0+372258 


120 


0.03 X 0.03 


130 


7 


1.2 X 10* 






S22 


J20273 1.7+372334 


80 


0.13x0.03 


150 


8 


1.5 X 10' 






S23 


J202731. 8+372154 


110 


0.24 X 0.03 


310 


16 


1.2 X 10' 






S24 


J202734.7+372311 


140 


0.18x0.09 


420 


22 


4.9 X 10"* 






S25 


J202923.6+401110 


120 


0.05 X 0.03 


140 


7 


6.1 X 10' 






S26 


J202924.8+401118 


880 


0.14x0.10 


2340 


<80-124 


<2.1-3.4xl0' 


1023.40 


AFGL 2591, cm'' free-free, 

maser'* (H2O/OH?) 


S27 


J20293 1.8+3901 13 


130 


0.19x0.09 


420 


< 14-22 


<2.7-4.2xl0'* 


23.55 


MSX 77.9550+0.0058, 
cm' free-free 


S28 


J202936.8+390112 


120 


0.12x0.08 


270 


< 9 - 14 


< 3.6- 5.6 X 10"* 


204.20 


IRAS 20277+3851, cm' free-free 


S29 


J202958.3+401558 


160 


0.23 X 0.12 


710 


37 


3.4 X lO'' 






S30 


J203 112.6+400316 


360 


0.16x0.16 


1560 


82 


7.7 X 10" 




close to IRAS 20293+3952 


S31 


J203 114.2+400305 


110 


0.13x0.07 


220 


12 


6.5 X 10* 






S32 


J203120.3+385716 


250 


0.21 X 0.11 


1010 


54 


6.0 X 10* 






S33 


J203144.5+385639 


80 


0.14x0.11 


230 


12 


2.7 X 10* 






S34 


J203157.8+401830 


150 


0.20x0.12 


610 


32 


3.4 X 10* 




in IRDC G79.3-P3 


S35 


J203158.8+385836 


200 


0.12x0.06 


410 


22 


1.4 X 10' 


5.00 




S36 


J203220.8+401949 


80 


0.07 X 0.03 


100 


5 


2.5 X 10' 




in IRDC G79.34+0.33 


S37 


J203222.0+402012 


240 


0.18x0.12 


860 


45 


6.4 X 10* 


4.44 


IRAS 20305+4010 
in IRDC G79.34+0.33 


S38 


J203222.4+401919 


90 


0.23 X 0.12 


410 


22 


2.0 X 10* 




in IRDC G79.34+0.33 


S39 


J203229.0+401600 


170 


0.14x0.09 


440 


<15-23 


<4.2-6.6xl0* 


37.71 


MSX 79.2963+0.2835, 

cm' free-free 


S40 


J203232.2+401625 


80 


0.16x0.10 


240 


<8-13 


< 1.6-2.6x10* 


73.03 


MSX 79.3070+0.2768, cm' free-free 
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Table 1. Continued 





Fragment 
name 


g peak ^ 
1.2 mm 
(mJy/ 
beam) 


FWHM" 
(pc X pc) 


g int a 
1.2 mm 
(mJy) 


^1.2 mm'' 
(Mo) 


< >' 
(cm-3) 


MSX-i 
21 ytim 

(Jy) 


Name, comments, and coincidence 
with stellar activity signatures 


S41 


J203233.9+401653 


220 


0.19x0.11 


800 


27-42 


3.7-5.9x10^ 


49.4 


part of IRAS 20306+4005 


S42 


J203239.9+384603 


90 


0.18x0.13 


340 


18 


2.0 X lO'' 






S43 


J203240.8+384631 


140 


0.23 X 0.13 


650 


35 


2.9 X lO"* 


Abs 




CygX-NW region 


NWl 


J201938.8+405638 


280 


0.17x0.15 


1110 


<38-59 


<4.2-6.6xl0* 


400.66 


IRAS 20178+4046, cm' free-free 


NW2 


J201940.5+405704 


200 


0.23 X 0.07 


680 


36 


7.7 X 10" 






NW3 


J202029.9+412144 


80 


0.15 X 0.07 


190 


10 


4.2 X 10" 






NW4 


J202030.0+4 12207 


130 


0.07 X 0.05 


180 


10 


2.2 X 10' 






NWS 


J202030.7+412125 


420 


0.07 X 0.06 


640 


22-34 


3.2-5.0x10' 


85.85 


EM* LkHA 225S, maser'" (H2O) 


NW6 


J202032.1+412354 


130 


0.15 X 0.03 


250 


13 


2.0 X 10' 






NW7 


J202032.2+412123 


80 


0.07 X 0.03 


100 


6 


2.5 X 10' 






NWS 


J202323.0+411749 


90 


0.10x0.04 


150 


8 


1.2 X 10' 






NW9 


J202324.0+411739 


230 


0.10x0.04 


370 


19 


3.0 X 10' 


27.47 


IRAS 20216+4107 


NWIO 


J202343.6+411657 


160 


0.20 X 0.08 


520 


27 


5.7 X 10" 






NWll 


J202413.0+411700 


110 


0.13 X 0.09 


290 


15 


4.7 X 10" 






NW12 


J202414.3+421143 


140 


0.06 X 0.06 


200 


11 


2.1 X 10' 


Abs 




NW13 


J202419.5+421545 


100 


0.17x0.09 


290 


15 


3.5 X 10" 






NW14 


J20243 1.7+420423 


420 


0.12x0.09 


940 


32-50 


1.2-2.0x10' 


15.58 


Moll 24 = IRAS 20227+4154, 



maser" (H2O) 

Notes: 

(a) Peak flux, deconvolved FWHM size, and integrated flux derived from a 2D-Gaussian fit to the 1 .2 mm map after background 
subtraction (cf. Sect. l3.2l and Appendix B). An upper limit FWHM of 0.03 pc was assumed for unresolved fragments. 

(b) Mass derived from Col. 5 using Eq. ([TJ and assuming /fi.imm - 0-01 cm^g"' and T^ust = 20 K. A second estimate with 
r^just = 40 K is made for H 11 regions and the brightest infrared protostars of Cygnus X. 

(c) Volume-averaged density derived from Col. 4 and Col. 6 using Eq. 

(d) The 21 ;um fluxes are those given in the MSX point source catalog (MSX C6), except for CygX-N30, N43, N46, N47, N65 and 
S41, for which an aperture flux is measured on the MSX map. The MSX flux of two extended H n regions is equally distributed 
between CygX-N46 and N47, and between CygX-S18, S19, S20. All fragments detected at 21 jjm are also detected at 8 jum 
with the exception of CygX-N30, N46, N47, and N65. "Abs" means that sources are seen in absorption at 8 fim. 

(e) Observations of high-density tracers toward CygX-N4 suggest that this 1.2 mm source does not correspond to a dense cloud 
fragment. Its nature is yet unclear but it is probably not a site for future or ongoing star formation. 

Re ferences: (l )IKurtzet al. 1994i (2)IPestalozzi et al. 20051 (3)I Haschick et al. 198H (4)IMinier et al. 20011 (5) [Argon et al. 2000 
(6)|HarrisJJ7 3 ; (8) Braz & Epchtein 1983 ; (9) Zoonematkermani etal. 19901 (lO) IMolinari et al. 1998bt (1 D IPalla et al. 199 H (12) 
Pipheretal. 1976, (13) .Trinidad et al. 2003i (14) ,Pallaet al. 19951 
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Table 2. Properties of clumps (large-scale cloud structures) detected in the Cygnus X complex 





Structure 


FWHM^ 


5,'"t = 
1.2 mm 




<nH2 




Comments 




name 


(pc X pc) 


(mJy) 


(Mo) 


(cm" 






CygX-North region 


Cl-Nl 


J203534. 1+422004 


0.48 X 0.31 


8470 


665 


5.0 X 


10* 


contains CygX-Nl, N2, N3 


C1-N2 


J203605.0+420851 


1.46x0.94 


15150 


1189 


3.1 X 


10' 


Starless? 


C1-N3 


J203607.3+414001 


0.64 X 0.38 


7010 


550 


1.9 X 


10" 


contains CygX-N5, N6 


CI-N4 


J203651.5+413625 


0.66 X 0.65 


14780 


1161 


1.7 X 


10" 


contains CygX-NlO 


CI-N5 


J203657.9+421132 


0.58 X 0.49 


11690 


917 


2.5 X 


10" 


contains CygX-N12, N13 


C1-N6 


J20370 1.9+4 13453 


0.81 X 0.64 


13580 


1066 


1.2 X 


10" 


contains CygX-N14 


CI-N7 


J20371 1.9+413338 


1.45 X 0.90 


16370 


1285 


3.6 X 


103 


contains CygX-N15 


CI-N8 


J203728.4+421630 


1.80 X 1.01 


17060 


1339 


2.3 X 


103 




C1-N9 


J203729.8+421409 


1.32x0.73 


14350 


1127 


4.9 X 


103 


contains CygX-N18 


Cl-NlO 


J203805. 3+423937 


0.87 X 0.51 


10000 


785 


1.1 X 


10" 


Starless? contains CygX-N21, N22, N24 


Cl-Nll 


J203833.2+420619 


1.25 X 0.74 


14880 


1168 


5.5 X 


103 


Starless? contains CygX-N29 


Cl-Nl 2 


J203836.3+423937 


1.30x0.56 


10240 


804 


5.5 X 


103 


contains CygX-N28 


Cl-Nl 3 


J203836.8+423734 


0.48 X 0.40 


46930 


3684 


1.9 X 


10^ 


contains CygX-N30, N31, N32 


C1-N14 


J203900.5+422238 


0.88 X 0.44 


92460 


7259 


1.3 X 


10^ 


contains CygX-N36, N38, N41, N44, N48 


Cl-Nl 5 


J203901.3+421935 


0.56 X 0.39 


70210 


5513 


2.3 X 


10^ 


contains CygX-N42, N46, N47 


Cl-Nl 6 


J20390 1.7+422507 


0.91 X 0.60 


47360 


3718 


3.8 X 


10" 


contains CygX-N37, N43, N51, N53, N54 


Cl-Nl 7 


J203927.2+4 12009 


0.78x0.59 


15980 


1254 


1.7 X 


10" 


contains CygX-N58, N59 


Cl-Nl 8 


J203936.4+4 12001 


0.91 X 0.36 


5710 


449 


1.0 X 


10" 


Starless? contains CygX-N60 


Cl-Nl 9 


J203955.4+420020 


1.79 X 1.10 


18110 


1422 


2.2 X 


103 


contains CygX-N61 


C1-N20 


J204003.6+412814 


1.03x0.77 


14080 


1105 


6.6 X 


103 


contains CygX-N62 


C1-N21 


J204005.1+413210 


0.31 X 0.20 


3350 


263 


7.0 X 


10" 


contains CygX-N63 


C1-N22 


J204027.8+415707 


0.68 X 0.37 


5200 


408 


1.4 X 


10" 


contains CygX-N64, N65 


C1-N23 


J204033. 3+415045 


0.80 X 0.57 


8700 


683 


9.4 X 


103 


Starless? contains CygX-N67, N69 


CygX-South region 


Cl-Sl 


J202038.5+393753 


0.51 X 0.34 


13920 


1093 


6.3 X 


10" 


contains CygX-S6, S7, S8, S9 


C1-S2 


J202044.8+393525 


0.67 X 0.45 


4870 


382 


9.4 X 


103 


contains CygX-SlO 


C1-S3 


J202716.1+372231 


1.23x0.69 


26640 


2092 


1.1 X 


10" 


contains CygX-S15, 816, 817 


C1-S4 


J202726.9+372248 


0.73 X 0.53 


33420 


2624 


4.6 X 


10" 


contains CygX-818, 819, 820, S21 


C1-S5 


J202924.9+401116 


0.42 X 0.38 


10090 


792 


5.2 X 


10" 


contains CygX-825, 826 


C1-S6 


J202959.4+401555 


1.25 X 0.60 


11560 


908 


5.9 X 


103 


Starless? contains CygX-S29 


C1-S7 


J203 112. 1+4003 14 


0.62 X 0.49 


10710 


841 


2.1 X 


10" 


contains CygX-830, 831 


C1-S8 


J203 119.2+385731 


0.86 X 0.56 


8740 


686 


8.6 X 


103 


contains CygX-832 


C1-S9 


J203156.7+401833 


1.39x0.60 


12140 


953 


5.2 X 


103 


Starless? IRDC G79.33P3, contains CygX-834 


Cl-SlO 


J203221.9+401955 


0.92 X 0.79 


17490 


1373 


9.3 X 


103 


contains CygX-836, 837, 838 


Cl-Sll 


J203232.3+401635 


0.99 X 0.67 


14010 


1100 


8.5 X 


103 


contains CygX-839, 840, 841 


C1-S12 


J203240.7+384614 


0.67 X 0.52 


6480 


509 


1.0 X 


10" 


Starless? contains CygX-S42, S43 


CygX-NW region 


Cl-NWl 


J201939.0+405651 


0.73 X 0.64 


14190 


1114 


1.5 X 


10" 


contains CygX-NWl, NW2 


C1-NW2 


J202030.5+412132 


0.64x0.41 


6890 


541 


1.7 X 


10" 


contains CygX-NW3, NW4, NW5, NW7 


C1-NW3 


J202344.5+411701 


0.97 X 0.71 


11480 


901 


6.5 X 


103 


Starless?, contains CygX-NWlO 


C1-NW4 


J202420.3+421527 


0.96 X 0.66 


8160 


640 


5.3 X 


103 


Starless?, contains CygX-NW13 


C1-NW5 


J20243 1.7+420419 


0.73 X 0.60 


8670 


681 


9.8 X 


103 


contains CygX-NW14 



Notes: 



(a) Deconvolved FWHM size and integrated flux derived from a 2D-Gaussian fit to the 1 .2 mm map smoothed to a 55" beam. 

(b) Mass derived from Col. 4 using Eq. ([TJ and assuming /cj 2 mm - 0-01 cm^ g"' and Tjust = 15 K. 

(c) Volume-averaged density derived from Col. 3 and Col. 5 using Eq. (|2]l. 
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Fig. 12. MAMBO maps of CygX-North (left: gray-scale and contours, right: contours overlaid on 8 /im images obtained by MSX 
and converted to Jy sr ') extracted from Fig.|2^. Regions shown are south and north of DR17 (a-b) and around DR20 (c). The 
1.2 mm and 8 fim images have 11" and 20" angular resolutions, respectively. The compact cloud fragments discovered in MAMBO 
images (see Table [TJ are labeled and marked by crosses in the gray-scale plot. The infrared sources which coincide with a MAMBO 
cloud fragment are also indicated. Contour levels are logarithmic and go from 40 to 800 mJybeam"' in a and b, and from 75 to 
800 mJy beam ' in c. 



28 



F. Motte et al.: High-mass star formation in the Cygnus X complex 




Fig. 12. (continued) Same caption for the region around W75N (d) and DR23 (e), and south of DR23 (f). Contour levels are 
logarithmic and go from 40 to 4 800 mJy beam"' in d, and from 40 to 800 mJy beam"' in e and f. 
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Fig. 13. Same convention as Fig. [12] for MAMBO maps of CygX-South extracted from Fig. Regions shown are the most 
western part of CygX-South (a), around Moll21 (b), around AFGL 2591 (c), and south of DR5 (d). Contour levels are logarithmic 
and go from 60 to 180 mJybeam"' in a, from 60 to 2 000 mJybeam"' in b, from 40 to 200 mJybeam"' in c, and from 40 to 
1 000 mJy beam ' in d. 
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Fig. 13. (continued) Same caption for the region around DR9 (e) and DR13 (f). Contour levels are logarithmic and go from 40 to 
200 mJy beam"' in e, and from 60 to 300 mJy beam"' in f. 
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S -CygX-N¥ 




Fig. 14. Same convention as Fig. [T2]for MAMBO maps of CygX-NW extracted from Fig. |2j;. Regions shown are those of Fig. |2j;4 
(a), of Fig. |2};3 (b), and of Fig. |2};2 (c). Contour levels are logarithmic and go from 40 to 300 mJybeam"' in a, from 60 to 
400 mJy beam ' in b, and from 60 to 250 mJy beam ' in c. 
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Fig. 14. (continued) Same caption for the region around Moll24 (d). Contour levels are logarithmic and go from 60 to 
180 mJybeam ' in d. 
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Fig. 15. SiO(2 - 1) emission toward 12 additional fragments, ordered by decreasing Mi. 2mm- 



